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ABSTRACT 

We present new high angular resolution near-infrared spectroscopic observations of the nuclear 
star cluster surrounding the Milky Way's central supermassive black hole. Using the integral-field 
spectrograph OSIRIS on Keck II behind the laser-guide-star adaptive optics system, this spectroscopic 
survey enables us to separate early-type (young, 4-6 Myr) and late-type (old, > 1 Gyr) stars with a 
completeness of 50% down to K' = 15.5 mag, which corresponds to ~ 10 M Q for the early-type stars. 
This work increases the radial extent of reported OSIRIS/ Keck measurements by more than a factor 
of 3 from 4" to 14" (0.16 pc to 0.56 pc), along the projected disk of young stars. For our analysis, we 
implement a new method of completeness correction using a combination of star-planting simulations 
and Bayesian inference. We assign probabilities for the spectral type of every source detected in deep 
imaging down to K' = 15.5 mag using information from spectra, simulations, number counts, and the 
distribution of stars. The inferred radial surface-density profiles, S(i?) oc i?~ r , for the young stars 
and late-type giants are consistent with earlier results ( T ear i y = 0.93 ±0.09, Ti ate = 0.16 ±0.07). The 
late-type surface-density profile is approximately flat out to the edge of the survey. While the late-type 
stellar luminosity function is consistent with the Galactic bulge, the completeness-corrected luminosity 
function of the early-type stars has significantly more young stars at faint magnitudes compared to 
previous surveys with similar depth. This luminosity function indicates that the corresponding mass 
function of the young stars is likely less top-heavy than that inferred from previous surveys. 

Subject headings: Galaxy: center, infrared: stars, stars: early-type, stars: stars: luminosity function, 
techniques: high angular resolution, techniques: spectroscopic 



1. INTRODUCTION 

The star cluster at the center of the Milky Way has 
been observed extensively in the past owing to its unique 
position in the closest galactic nucleus. The study of its 
properties has led us to unique insights about its stel- 
lar population and has demonstrated the existence of 
a supermassive black hole at the Gala c tic center (e.g. 
Paumard et all [2001 iGhez et all [200l iGillessen et all 
2009f l. The nuclear star cluster is composed mainly of a 
ma ssive old stella r cluster with a half-light radius of 5-10 
pc (|Schodcl 20TTJ) . At the center of the cluster, located 
within the central ~ 0.5 pc, is a concentration of young 
stars (of age 4-6 Myr old) that dominates the luminosity 
of this region. These two components provide us with 
different probes of the physical conditions near a super- 
massive black hole. The presence of the young stars in 
the strong tidal field of the black hole allows us to study 
star formation in an extreme environment and provides a 
test of the universality of the initial mass function. The 
late-type old stars, on the other hand, provide us with 
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a test of the long-term interactions between a star clus- 
ter and a supermassive black hole. This has implications 
for black hole growth as well as the inward migration of 
compact objects. 

Our understanding of the nuclear star cluster in the 
Galactic center is driven in large part by progressively 
more advanced observing capabilities. Seeing-limited 
observations of the Galactic center in the infrared en- 
abled the identification of the nuclear star cluster as 
a peak in th e stellar density toward the center of the 
Galaxy (e.g. iBecklin fc Neugebauerl Il968f ). Subsequent 
spectroscopy led to the discovery that the center of the 
cluster also hosts a n umber of bright emission line stars 
l|Krabbe et al.l Il991f) . indicating that a population of 
young stars resides within the central ~ 0.5 pc. However, 
because of the high density of stars in this region, it was 
not possible to disentangle the two populations of stars 
through seeing limited observations. Spectroscopy of the 
spatially integrated light showed a decrease in CO equiv- 
alent width toward the center of the cluster, which can 
be due to either a decrease in the number of red giants 
or contamination of the sp ectra by the bright Wolf-Rayet 
(WR) stars in the region (jHaller et al.lll996| ). 

These limitations were greatly alleviated by the ad- 
vent of adaptive optics, which allowed diffraction-limited 
imaging and spectroscopy in the near-infrared on 8-10 m 
class telescopes. Adaptive optics imaging enabled mea- 
surements of the number counts of stars as well as their 
proper motions in the p lane of the sky (e.g. IGenzel et al.1 
l2003tlGhez et al1l20"0 5). Integral-field spectroscopy pro- 
vided the crucial ability to separate the population of 
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young stars from that of the old red giants, thus en- 
abli ng the study of th e two populatio ns independently 
(e.g. lEisenhauer et al.ll2005t IDo et alj|2009al) . In terms 
of studying the characteristics of star and cluster forma- 
tion in this region, these advances provide two key ob- 
scrvables: the surface-density profile and the luminosity 
function of the cluster. 

The surface-density profile provides one of the ob- 
servable features of the dynamical state of the cluster. 
Early in the formation of the cluster, the stellar dis- 
tribution reflects its origin; for example, about half of 
the young stars are observed to be distributed in a thin, 
clockwise-rotating stellar disk with a steep projected ra- 
dial surface-density p r ofile of ~ 1/R 2 in the disk plane 
(jPaumard et all 120061; ILu et alj[2009h . This may be in- 
dicative of their in-situ form ation in an accretion disk 
(jLevin fe B cloboroclov 2003). On the other hand, over 
time, the cluster will become dynamically relaxed with 
respect to the black hole and settle into a steady-state 
density profile, with all traces of its origin removed. 
iBahcall fc Wolj (|1977|) predicted that star clusters with 
a massive black hole should contain a cusp with a spatial 
density profile of r~ 7 / 4 to r~ 3 / 2 , depending on whether 
the cluster has a single mass population or contains mul- 
tiple mass components. This property helped to fa- 
cilitate calculations such as the growth of black holes 
by stars and the in-spiral rate of compact stellar rem- 
nants in galactic nuclei, as the power-law exponent is 
one the most uncertain parameters for describing th e 
distribution of stars (e.g. iPreto fc Amaro-Seoandl20l"ol ). 
While the red giants at the Galactic center (> 1 Gyr) 
may have had time to dynamically relax, they unexpect- 
edly show a core-like (i.e., flat) surfa ce-density profile 
(|Buchholz et alj|2009t IDo et al.ll2009af) . At present, the 
origin of the flattening of the surface-density profile is 
unclear. The flat core could arise from secular evolution 
of the cluster, such as resonant relaxation, or caused by a 
drast ic event such as the in-fall of another massive blac k 
hole (|Merrittll20ldt IMurphvll201lt IMadiganltalll20Tl . 
Increasingly refined measurements of the properties of 
the structure of the old stellar population, such as its 
spatial density profile and core radius, are necessary to 
make progress. 

The luminosity function is one of the most fundamen- 
tal observable parameters of any stellar population. It 
is a measure of the relative distribution of stellar lu- 
minosities and can be used to determine such proper- 
ties as the age, star formation history, and initial mass 
function (IMF) of the cluster. Much of the early work 
on the near-infrared luminosity function at the Galac- 
tic center was aimed at understanding the old popula- 
tion of bright giants, as many of them can be spatially 
resolved with seei ng-limited imaging and spectroscopy. 
IBlum et al.l (|1996l ) conducted one of the the most com- 
plete near-IR photometric surveys of the central 2' (~ 5 
pc) of the Galactic center possible under seeing-limited 
conditions. Their observed luminosity function reached 
K ~ 12.5, where stellar crowding started to dominate. 
They found that down to these magnitudes, the K lu- 
minosity function in t his re gion is consistent with that 
found bv lTiede et al.l (|1995i ) for Baade's Window, a low 
extinction region several degrees from the Galactic cen- 
ter. While this suggests that the Galactic center may 



have the same star formation history and composition as 
the inner bulge, the observations were not deep enough 
to reach the red clump at K = 15.5, where most of the 
red giants at the Galactic center are manifested. 

The IMF is one of the most important observational 
signatures that connect star formati on theories with ob- 
servations (|McKee fc Ostrikerl 120071 ). As most observa- 
tions from the local universe show a remarkably consis- 
tent IMF across different star formation environments, 
there is substantial interest in whether the stellar IMF 
is universal, especially in extreme environments like th e 
Galactic center (see review from iBastian et al.1 12010T ) . 
The best population for constraining the IMF of stars at 
the Galactic center lies in the young stars within the cen- 
tral parsec, which have recently become observationally 
acce ssible through integral-field spectroscopy behind AO 
(e.g. lEisenhauer et alJl2005b IDo et alJl2009ah. Us i ng the 
AO-fed IFU SINFONI at VLT. iPaumard et al.l ([20061 ) 
were the first to construct a K luminosity function from 
a spectroscopically selected sample of young stars. Their 
sample was largely limited to stars brighter than K = 13 
mag, which is near the transition between evolved OB 
stars and main-sequence B stars. This corresponds to 
measuring the mass function only for the evolved mas- 
sive stars. It is necessary to observe the main sequence 
for reliable mass function measurements as the stellar 
atmosphere and evolutionary models may have large un- 
certa in ties for the massi ve evolved stars (jMartins et al.1 
l2007h . IDo et al.l (|2009af ) went deeper with the OSIRIS 
spectrograph on Keck II to K < 15.5, which provided a 
sample of the early-type main sequence B stars. How- 
ever, they did not attempt to derive the mass function 
from the observed luminosity function. iBartko et al.l 
(|2010( ) also achieved a similar depth for spectroscopic 
observations in fields sampling out to ~ 1 pc from the 
center, largely perpendicular to the disk of young stars. 
They used star planting simulations to derive a com- 
pleteness correction for the K luminosity function within 
projected radius 0."8 < R < 12" (0.03 to 0.5 pc) and 
through stellar population synthesis modeling, concluded 
that the young stars have a very top-heavy IMF, with 
dN /dm oc to~ 45± 0,3 , compared to a Salpeter IMF of 
dN I dm cc m " 2 - 35 (jBartko et alJl20Tol ). In comparison, 
the young stars in the central 0."8 (sometimes called the 
S-stars) show a slope that is consistent with Salpeter, 
with T = -2.15±0.3 for stars with K > 14.0 mag. There 
are suggestions that the S-stars may not originate from 
the same star formation event that formed the young 
stars further due to the lack of stars more massive than 
ea rly B main sequen ce stars in this region (see review 
in lGenzel et al.ll2010T l . Some theoretical studies however 
suggest that the S-stars might have originated further 
out, and wer e brought in by a com bination of dynamical 
events fe.g.. iMadigan et al.ll2009T ). One challenge in in- 
terpreting these luminosity functions is that at a depth 
of K' < 15.5, there is only a limited magnitude range 
from the start of the main sequence at K' ss 14.0; this 
limited range means that the faintest magnitude bin has 
a large impact on the slope of the luminosity function, 
and hence the mass function. 

Because of the importance of the faint end of the lu- 
minosity function, it is crucial to understand and care- 
fully account for the assumptions inherent in complete- 
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ness correction. Completeness correction attempts to 
characterize an underlying population, in which only a 
certain number of sources can be observed or identified. 
In the case of the luminosity function of young stars at 
the Galactic center, the aim would be to quantify the 
number of young stars, given that not all sources will 
have spectral identification. One s tandard approach is 
through star-planting simulations as lBartko et al.l (|2010f ) 
have done in order to characterize the recovery rate of 
stars at a given magnitude; the completeness-corrected 
count is then derived by dividing the observed number of 
stars by the fraction of undetected sources. This method 
utilizes no other information than that provided by the 
star counting simulation. However, for the case of the 
Galactic center, we have much more information on the 
underlying population that can be used in the complete- 
ness correction: (I) the sources in the magnitude range 
where spectroscopy is possible have been extensively im- 
aged and their counts are nearly 100% complete; this 
means that at a minimum, the completeness-corrected 
number of sources should not exceed the number of ex- 
isting sources. (2) Given the luminosity function of stars 
at the Galactic center at the current spectroscopic sen- 
sitivity, we are mainly sensitive to two types of sources 
- old late-type giants and young stars (Figure [T]). (3) 
The surface density profiles of the early and late-type 
stars also give information on their expected fraction as 
a function of projected distance from Sgr A*. 

Here, we present new integral-field spectroscopic ob- 
servations that extend to a distance of ~ 0.5 pc from Sgr 
A* in the direction along the projected major axis of the 
plane of the young stella r disk, with a sky position an- 
gle (P A) of 105 deg (e.g. iPaumard et"al|[200l iLu et al.l 
|2009( ). Previously, spectroscopic coverage beyond about 
0.25 pc in this region was l imited to lower spati al res- 
olution measurements (e.g. IPaumard et all l2006f ). We 
obtain spectra for about 400 stars with K' < 15.5 mag, 
which allow us to investigate the radial profile and lu- 
minosity functions of both the early-type (young) and 
late-type (old) populations. In Sections [2] and El we de- 
scribe the new observations and data reduction, while in 
Section |4] we review our method for assigning spectral 
types. In Section |4~31 we adapt a method from Bayesian 
inference with star-planting simulations to infer proba- 
bilities for the spectral types of all sources with K' < 15.5 
and to establish the spectral completeness of our survey. 
This incorporates all available information from spectra, 
number counts, and knowledge of the radial distribution 
of stars at the Galactic center to estimate the stellar 
population. In Section [51 we present the resulting sur- 
face number density profiles and K' luminosity functions 
of the early and late- type stars and in Section |6l we dis- 
cuss the implications of these results for the mass func- 
tion of the young stars and their origins, and we present 
the implications for the radial d istribution of o ld stars. 
In a companion paper (Paper II. lLu et al.l[2012h we will 
derive the mass function of the young stars. Section [7] 
summarizes our conclusions. 

2. OBSERVATIONS 

Observations of the central 0.5 parsec of the Galaxy 
consist of (1) spectroscopy to distinguish young, hot stars 
from old, cool giants and (2) photometry to measure the 
brightness of each young star. Both the spectroscopic 
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Figure 1. Theoretical color- magnitude diagram of the Galac- 
tic center in the NIR showing the expected observed K' magni- 
tude for stars of different spectral types behind 2.7 magnitudes 
of extinction at K' and at 8 kpc (the H-K colors are intrinsic 
colors). The locations at 50% completeness for spectroscopy and 
imaging are also shown. Most of the stars observable with spec- 
troscopy in this region are either young stars (blue) or late-type 
giants (re d). The colors and magnitudes in this plot are de - 
rived fr om Williams & Antonopoulou ( 19811): iDucati et all l|2001l '): 
Martins k, Plcz (2006); Wcencr (2007)1 

and imaging observations were obtained in conjunction 
with the laser-guide-star adapti ve optics (LGS AO) sys- 
tem on the Keck II telescope (|Wizinowich et ahl 120061 : 
Ivan Dam et al.l 120061 ): the laser guide star was propa- 
gated at the center of the field of view for each observa- 
tion, and for low-order tip-tilt corrections, we used the 
R=13.7 mag star, USNO 0600-28577051, which is located 
~19" from SgrA*. Details specific to the spectroscopic 
and imaging observations are described in Sections 12.11 
and 12.21 respectively. 

2.1. Spectroscopy 

Near-IR integral-field spectra of the Galactic center 
were obtained between 2007 and 2011 using the OH- 
Suppressing Infra- Red Imaging Spectrograph (OSIRIS, 
ILarkin et al.l I2006h . The primary observations for this 
work constitute a survey through the narrow-band fil- 
ter Kn3 (2.121 to 2.220 urn) This in c ludes b oth ob- 
servations initially reported in IDo et al.l (|2009al) and iDol 
(120101) and new 2010-2011 observations, which increase 
the radial extent of this survey by more than a factor 
of 3 from 4" to 14" (0.16 pc to 0.56 pc, see Figure 
[2]). We refer to the combination of the original survey 
and this new extension of our survey as the Galactic 
Center OSIRIS Wide-field Survey (GCOWS). Our ini- 
tial work, which covered a 8" x 6" region centered on 
Sgr A*, used OSIRIS's 35 mas plate scale (field of view 
of 1.58" x 2.24"). The new GCOWS observations are lo- 
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cated at larger projected distances from the Galactic cen- 
ter than the previous work and were obtained with a 50 
mas plate scale (field of view of 2.25" x 3.2"); as the stel- 
lar densities are lower in this region, this provides a good 
compromise between spatial resolution and field of view. 
The new GCOWS fields cover a region of approx imately 
10" x 7.2" east of the survey reported inlDol (|2010l) . along 
the major axis of the projected disk plane of the clockwise 
dis k of young star s at a PA of 105 degrees, as mea sured 
bvlLu et all (12009ft (s e e alsolLevin fc Beloborodoyll2003l : 
iPaumard et all 120061 : iBartko et all 120091 : lYeldal 120121 ). 
Each of the new fields is observed with a six-point dither 
pattern of 900 s per frame, in which the dithers have 
small (~0".l) offsets from one another. The larger plate 
scale allows us to reach sensitivity comparable to that of 
the previous 35 mas plate-scale observations, which have 
about 9 dithers per field. The total s urface area of al l 
the observations (including those from lDo et alJl2~009a[ ) 
is 113.7 square arcseconds. We also observe 7 pointings 
(of various total integration times) within the Kn3 survey 
region with the K broad-band filter (Kbb, 1.965-2.382 
^m) in either the 35 mas or 50 mas plate scale, depend- 
ing on the stellar density. The broad-band observations 
are used to verify the spectral-types of a sample of stars 
see Appendix ©. Table [U summarizes the details of the 
complete survey, including: field locations, integration 
times, dates of observations, and data quality. 

For calibration purposes, we observe skies after the 
Galactic center observations. These observations arc 
used to determine the stability of the wavelength solution 
with the OH sky lines. Sky subtractions for the Galactic 
center spectra are done using local sky measurements in 
each of the science data cubes as described in Section [3l 
To remove atmospheric telluric absorption lines, we also 
observe an A0V (HD195500 or HD155379) and a G2V 
(HD193193 or HD150437) star each night. 

2.2. Imaging 

Photometric observations for the individual stars in the 
central parsec of the Galaxy were conducted using the K'- 
band (K') filter at A = 2.12 /im (AA = 0.35 (im) in order 
to identify stars, measure their positions and K' bright- 
ness, and estimate the completeness of our star counts. 
Observations were taken in 2006, 2008, and 2010 with 
the NIRC2 instrument (PI: K. Matthews). The NIRC2 
field of view is ~10 " with a pixel scale of 9.95 mas/pixel 
(|Yelda et alJl2010h . To image all of the young stars in 
this region, a mosaic was constructed covering 27" x 27" 
roughly centered on Sgr A*. Individual exposures at each 
pointing had an integration time of ti nt = 28 s (2.8 s x 
10 coadds). The mosaic dither pattern always consisted 
of a 3x3 position box pattern with 8"5 steps between 
each position and with ~ 4 exposures used at each po- 
sition^. For 2 of the 3 epochs, we also completed a 2x2 
position box pattern with 4" steps, and 2 exposures at 
each position to provide large overlaps between all the 
tiles in the mosaic. 

3. DATA REDUCTION 

3.1. Spectroscopy 

8 The observations were typically done under less than ideal 
seeing conditions and some individual exposures were rejected due 
to poor AO correction. 



Data reduction and extraction of ne w spectra were per - 
formed in a manner similar to that of lDo et"ail (|2009aft . 
For these data, we used version 2.3 of the OSIRIS 
pipeline, as provided by the instrument team. This ver- 
sion includes a new wavelength solution for the instru- 
ment in 2009, which was subject to changes in tempera- 
ture at that time. This wavelength solution was verified 
by comparing the locations of OH sky emission lines. 
The pipeline also removes electronic cross-talk, corrects 
cosmic rays, and assembles the data cubes. 

Stellar spectra are extracted from the GCOWS Kn3 
observations of all stars brighter than a differential ex- 
tinction corrected <15.5 mag within our field of 
view, as identified from deep NIRC2 imaging (Section 
[372]) . To register the OSIRIS data cubes to the LGS 
AO K 1 images, the point-spread functio n (PSF) fitting 
routine, StarFinder (jDiolaiti et al.l I2000T ). is run on an 
image produced from collapsing the OSIRIS data cube 
along the spectral dimension, producing an OSIRIS star 
list. Positions from StarFinder are matched to positions 
derived from LGS AO K 1 images in the same region. Be- 
cause the images have higher spatial resolution and bet- 
ter PSF characterization than the OSIRIS observations, 
they allow us to identify the locations of stars that may 
have been missed by StarFinder on the OSIRIS cubes. 
Spectra are then extracted with a circular aperture cen- 
tered at the location of each star (K^ < 15.5) detected 
in imaging at each spectral channel. We use an aperture 
radius between 1 and 2 pixels (50-75 mas), depending on 
the distance to the nearest source. For sky and back- 
ground subtraction, we use the median flux values in an 
annulus with an inner radius of 1-2 pixels and an outer 
radius of 2-4 pixels. To remove atmospheric telluric lines, 
we divide the spectra by that of a blackbody-removed 
spectrum of an A star each night. The A star is feature- 
less in the wavelength-region of interest except for the 
strong Br 7 line, which we replace by using the spectrum 
of a G star calibrator divided by the solar spectrum over 
the region 2.155-2.175 ^m. Within the entire GCOWS 
data set, we ext ract a tot al of 400 spectra, including 
those reported in lDol ()2010f ). which are re-extracted and 
re-analyzed here. 

Stellar spectra are extracted from the Kbb observa- 
tions for 12 stars (described in more detail in Appendix 
[A"f . For Kbb data taken in the 50 mas plate scale, stellar 
spectra are extracted in a similar manner to the extrac- 
tion described above. For the Kbb data taken in the 35 
mas plate scale, an aperture radius of 2 pixels is used, 
with a sky annulus defined from 2-4 pixels. 

3.2. Photometry 

Each tile of the NIRC2 photometric mosaic is reduced 
separately. This is necessary because the AO point 
spread function (PSF) varies with time and position and 
the correct PSF is required for precise photometry. Our 
NIRC2 data reduction pipeline is used to subtract dark 
current and sky emission, flatten the field, remove bad 
pixels and cosmic rays, and apply corrections for instru- 
ment al and atmospheric distortion fLu 2008; Yeld a et al.l 
120101 ). For each tile in the dither pattern, the individual 
exposures at that pointing were combined. Additionally, 
three subset- images are created for each tile with 1/3 of 
the exposures in order to estimate uncertainties. 

Stellar photometry and astrometry are extracted using 



■5 



Table 1 



Summary of OSIRIS observations 


Field Name 


Field Center 3 - 


Date 


N frames * tint 


Plate Scale 


FWHM b 


Filter 


Published 


PA 




(") 


(UT) 


00 


(mas) 


(mas) 






(degrees) 


GC Central (C) 


0, 


2008 May 16 


11 x 900 


35 


84 X 85 


Kn3 


1 


285 


GC East (E) 


2.88, -0.67 


2007 July 18 


10 x 900 


35 


85 x 70 


Kn3 


1 


285 


GC South (S) 


-0.69, -2.00 


2007 July 19 


10 x 900 


35 


73 x 63 


Kn3 


1 


285 


GC West (W) 


-2.70, 0.74 


2007 July 20 


11 x 900 


35 


110 x 86 


Kn3 


1 


285 


GC Southeast (SE) 


1.67, -2.23 


2008 June 03 


11 x 900 


35 


68 x 63 


Kn3 


1 


285 


GC North (N) 


0.33, 2.01 


2008 June 07 


7 x 900 


35 


102 x 85 


Kn3 


1 


285 




2008 June 10 


5 x 900 


35 


75 x 70 


Kn3 


1 


285 


GC Northeast (NE) 


2.55, 1.27 


2008 June 10 


5 x 900 


35 


74 x 68 


Kn3 


1 


285 


GC Southwest (SW) 


-2.9, -1.12 


2009 May 26 


4 x 900 


35 


92 x 80 


Kn3 


2 


285 


GC Northwest (NW) 


-1.99, 2.42 


2009 July 21 


6 x 900 


35 


71 x 64 


Kn3 


2 


285 


E2-1 


5.43, 0.99 


2010 May 6 


6 x 900 


50 


94 x 96 


Kn3 


3 


285 


E2-2 


4.8, -1.4 


2010 May 7 


6 x 900 


50 


88 x 79 


Kn3 


3 


285 


E2-3 


4.16, -3.75 


2010 July 28 


6 x 900 


50 


104 x 86 


Kn3 


3 


285 


E3-1 


8.59, 0.15 


2010 May 9 


6 x 900 


50 


79 x 86 


Kn3 


3 


285 


E3-2 


7.94, -2.21 


2010 May 7 


1 x 900 


50 


72 x 77 


Kn3 


3 


285 






2010 May 9 


5 x 900 


50 


79 x 86 


Kn3 


3 


285 


E3-3 


7.31, -4.57 


2010 July 29 


6 x 900 


50 


95 x 86 


Kn3 


3 


285 


E4-1 


11.73, -0.68 


2010 May 10 


6 x 900 


50 


97 x 84 


Kn3 


3 


285 


E4-2 


11.08, -3.04 


2010 May 9 


1 x 900 


50 


79 x 86 


Kn3 


3 


285 






2010 May 10 


5 x 900 


50 


97 x 94 


Kn3 


3 


285 


E4-3 


10.44, -5.41 


2010 July 29 


1 x 900 


50 


95 x 86 


Kn3 


3 


285 






2010 July 30 


5 x 900 


50 


106 x 90 


Kn3 


3 


285 


S2-1 


0.69, -4.16 


2010 August 1 


6 x 900 


50 


102 x 85 


Kn3 


3 


285 


Verification Field 3 


4.29, -2.20 


2011 July 25 


4 x 600 


50 


80 


Kbb 


3 


140 


Verification Field 4 


4.28, -3.76 


2011 July 25 


2 x 600 


50 


80 


Kbb 


3 


140 


Verification Field 5 


1.37, -2.38 


2011 August 17 


2 x 600 


50 


90 


Kbb 


3 


105 


Verification Field 6 


-1.34, -2.49 


2011 August 17 


3 x 600 


50 


90 


Kbb 


3 


105 


Verification Field 7 


5.52, 0.99 


2011 August 17 


1 x 600 


50 


90 


Kbb 


3 


105 


Verification Field 8 


-3.30, -0.65 


2011 August 17 


1 x 600 


50 


90 


Kbb 


3 


105 


Verification Field 9 


1.58, -0.62 


2011 August 25 


3 x 600 


35 


86 


Kbb 


3 


90 


Verification Field 11 


12.22, -4.00 


2011 August 25 


2 x 600 


35 


86 


Kbb 


3 


90 



a RA and DEC offset from Sgr A* (RA offset is positive to the east), 
k Average FWHM of a relatively isolated star for the night, found frc 
c 1 - IDo et alj I2009al) . 2 -[Dc] 1 12010) . 3 - this work 

StarF inder with the same setup described in lYelda et al.l 
(|2010ft . The resulting star lists are photometrically cali- 
brated using a samp le of stellar magnitudes reported in 
ISchodel et al.l (|2010f ). converted from the Ks filter to the 
K' filter as described in Appendix [Cj Uncertainties are 
estimated empirically for each tile by taking the error on 
the mean flux and position measurements from the tile's 
three subset-images. Sources not detected in the tile's 
combined image and 3 subset images are thrown out as 
spurious artifacts. Starlists for all the tiles are then mo- 
saicked together to create a single master starlist for each 
epoch covering the entire 27" x27". The photometry for 
stars in the mosaicked starlist is the error-weighted av- 
erage flux of all the tiles in which a star is present. The 
photometric errors are either the weighted standard de- 
viation of the fluxes in all the tiles, STD we i g hted(ft), or 
the average flux error, AVG(fj/ i ), whichever is larger. 

The mosaicked starlists from 2006, 2008, and 2010 are 
aligned together. Sources are dropped that are not de- 
tected in at least 2 of the 3 epochs, which throws out most 
spurious detections due to PSF artifacts and cosmic rays. 
Some stars may have intrinsic brightness variations, so 
we adopt, as our final photometric measurements, the 
time-averaged flux and RMS error, weighted by the flux 
errors at each epoch. Due to the small number of epochs 
used to estimate the error, we impose a minimum pho- 
tometric uncertainty of 0.02 mag. 

Photometry for stars at the Galactic center must also 
be corrected for strong and spatially variable extinction, 



a two-dimensional Gaussian fit to the source. 

even at near-infrared wavelengths. A detailed extinction 
map has been created for the regi on from near-infrare d 
photometry of red-clump stars by ISchodel et al.l (|2010|) . 
This extinction map is used to apply differential extinc- 
tion corrections to individual stars, thereby shifting all 
the observed stars in our NIRC2 imaging to a common 
extinction value of Ajf 5 = 2.7, the mean extinction value 
for the region. Before applying the differential extinction 
correction, we convert our observed K' magnitude to a 
Ks magnitude using filter conversions computed from a 
synthetic atmosphere with T e ff =30,000 K for early- type 
stars, and T e //=4,000 K for late-type stars and un-typed 
stars ( Appendix [Cj) . For the untyped sources, the error in 
assuming the wrong spectral type is less than the typical 
photometric error. After correcting for extinction, the 
Ks photometry is converted back to K' magnitudes and 
the differential extinction-corrected K^ photometry is 
used throughout the paper. 

The high stellar density and the large brightness con- 
trast of stars at the Galactic center cause some stars 
to be undetectable in the NIRC2 images. The imaging 
completeness as a function of position and brightness is 
estimated using star planting simulations described in 
Appendix IC.ll The average resulting completeness is 
94% at K' AA =15.5 and 41% at K' AA = 18 mag in the 
GCOWS field of view (these values are comparable for 
observed K'). 

4. SPECTRAL TYPING AND BAYESIAN INFERENCE 
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RA Offset (arcsec) 



Figure 2. The spatial coverage of the OSIRIS survey fields in the Kn3 filter (dashed blue) as well as the K broad-band spectral verification 
fields (dotted red). The survey is designed to increase the ra dial coverage along the orientation of the young stellar disk on the plane of the 
sky (dashed black, Q = 105 deg, i = 115 deg; ILu et aT|[2009l1 . The survey extends radially out to a projected distance of ~ 0.5 pc (dashed 
black) from Sgr A* (black diamond). The survey region from[Bartko et al. (2013, dot-dashed orange) overlaps the GCOWS fields closest 
to Sgr A*. 



For the purposes of this study, we wish to differentiate 
the Wolf-rayet (WR) and O /B stars (main sequence and 
evolved) from those of later spectral-types such as the 
evolved M and K giants. We will refer to early- type stars 
as stars with a spectral type of B or earlier (including 
the WR stars) and late-types as all stars with spectral- 
types later than B. As these two groups of stars were 
formed at very different times, we can use their early 
or late-type status as a proxy for age in measurements 
of the luminosity function and surface-density profiles of 
the two populations. In this section, we describe our 
method of classifying stars under the hypotheses that 
they are either early- type (He), or late- type (Hi). The 
goal is to assign each star a probability of being early- 
type, Pe, or late- type, Pl, and with the constraint that 
Pl + Pe = 1. This process is composed of the following 
steps: 

1. Manually assign each star as either early- type 
(Pe = 1), late-type (Pl = 1), or untyped using 
spectral classification criteria laid out in Section 

ED 

2. Use the sample of manually typed stars with K' 
> 14.0 to train the Bayesian algorithm to recognize 
the properties of early and late-type stars. This is 
accomplished by constructing the probability dis- 
tributions of Na I and Br 7 equivalent widths for 



both the early- type and late- type stars with K'>14 
(Section H2J). 

3. For all untyped sources, assign probabilities based 
on the Bayesian evidence for the early-type and 
late-type hypotheses using the above training sam- 
ple and extensive star planting simulations (Section 

on . 

For all stars we have extracted (K' AA <15.5), this analysis 
yields the probability that each star is either early-type 
or late-type. 

4.1. Manual Spectral Types 

We group the stellar spectra by eye into three groups: 
1) late-type, 2) early-type, or 3) untyped. Stars with 
significant Na I features are classified as late-type (219 
stars). The sources with Br 7 absorption and no Na I fea- 
tures are classified as early-type (44 stars) . Bright (K' < 
13.0) stars with featureless spectra between 2.121 to 
2.220 /im are also classified as early-type (23 stars); these 
sources are most likely O V or O/B I stars which can have 
very weak Br 7 absorption or emission (jHanson et al.l 
119960 . We also identify 12 WR stars, all of w hich were 
previously identified by iPaumard et all (|2006h . as early- 
type. The remaining stars with unclear spectral features 
are classified as untyped; all stars with SNR < 5 are also 
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Table 2 

Observed distribution of spectral line widths for 
K' > 14 



SpT 


Na EW 

(A) 


(A) 


Br 7 EW 

(A) 


(A) 


Late- type 
Early-type 


3.7 
-0.2 


1.3 
0.6 


1.9 
3.3 


1.2 
1.8 



classified as untyped. Th e above c r iteria a re slightly dif- 
feren t from those used in IDo et all (|2009a| ) . The revision 
from lDo et al.l (|2009a[ ) is based on the detection of a few 
yellow giants in the survey region, which have smaller 
Na I equivalent widths than the bul k of the K and M 
giants (some hav e been identified by iBlum et~all[200l 
IPfuhl et al.|[20TTI ). In the process of revising these crite- 
ria, we obtained Kbb spectra of a small subset of sources 
to verify that this spectral typing method is robust (see 
Appendix |A| . Figure [3] shows the locations of the 286 
stars with manually determined spectral classifications. 
Detailed properties for these late, non-WR early-type, 
and WR stars arc reported in Table HI Table [3 and Ta- 
ble [H respectively. 

The sensitivity of the fields observed in 2 010 with 
the 50 mas plate scale is similar to that from IDo et all 
(20093). Table [3] summarizes the completeness of each 
of the new fields compared to imaging in 0.5 mag bins. 
For the entire GCOWS sample, we are able to spectral- 
type about 50% of the sources known from imaging with 
K-aa between 15.0 to 15.5 mag (Figure [3]). The radial 
dependence of the spectroscopic completeness is shown 
in Figure [5] for K' = 14.5 — 15.5 mag. This is com pared 
to the imaging completeness presented in Appcndix lC.il 
The spectroscopic sensitivity drops dramatically around 
bright stars (Figure |3]). For example, no stars were 
spectral-typed at distances closer than 0."25 from the 
IRS sources. Key factors that contribute to the incom- 
pleteness of our observations are halo noise from bright 
stars, background gas emission lines, and crowding in the 
central regions. 

4.2. Training sample: observed line width distributions 

A key component of assigning probabilities to the un- 
typed stars is understanding the expected distribution of 
measured Na I and Br 7 equivalent widths of observable 
stars at the Galactic center. We use all manually typed 
stars with K' > 14.0 to construct a distribution function 
for the Na I and Br 7 equivalent widths of early-type 
and late- type stars (Figure 2]); these stars are chosen for 
the training sample because the majority of the untyped 
stars have K' > 14.0 mag. The method used to measure 
equivalent widths is described in Appendix |B] Typical 
errors in equivalent width of Na I for this sample for the 
early and late- type stars are about 0.5 and 1 A, respec- 
tively. The equivalent width errors for Br 7 area about 
0.5 and 0.7 A. A Gaussian was fit to each distribution of 
equivalent widths (the best fits are summarized in Table 
E}. These Gauss ian distributions are used as priors in 
the following Bayesian analysis of the untyped stars. 

9 As in IDo et al.l (20093) we exclude the star SO-32 from our 
analysis because it is a known foreground source. 



4.3. Statistical Spectral Types as an Approach to 
Completeness Correction 

Spectra of the untyped sources contain important in- 
formation about the relative completeness of early-type 
and late-type classifications, which has previously not 
been incorporated into analyses of luminosity functions 
and radial density profiles. In the following section, we 
develop a new approach to completeness correction that 
utilizes the spectra of untyped sources, along with ex- 
tensive star planting simulations, to assign each untyped 
source a probability of being either early-type or late- 
type. In essence, we compare two hypotheses: (1) the 
star is late-type (Hl) or (2) the star is early-type (He)- 
The goal is to compare the relative strengths of these two 
hypotheses and assign a probability for a given star to be 
early or late-type. To accomplish this goal, the Bayesian 
evidence is computed for both hypotheses given the ob- 
servations and our knowledge about the expected spec- 
tral features of these sources (Section l4.3.1[) . The relative 
strengths of the hypotheses is the ratio of the Bayesian 
evidence. To calibrate the Bayesian evidence and deter- 
mine probabilities, we perform extensive star planting 
simulations (Section l4.3.2[) and make use of two types of 
priors (Section I4.3.3[) . The results of these analyses are 
provided in Table [9] 

4.3.1. Bayesian Evidence 

First, we consider the Bayesian evidence, which is the 
likelihood of obtaining observed data, x, given a specific 
hypothesis, H , marginalized over all possible model pa- 
rameters, 9: 

P(x\H) = J P{x\9,H)P(9\H)d9 (1) 

P(x\9,H) is the likelihood of observing x for a specific 
model parameter 9, and P(9\H) incorporates prior infor- 
mation about the distribution of model parameters. In 
our case, the Bayesian evidence for the two hypotheses 
are 

P(Na obs ,Br obs \Na,Br,H E ) 

-OO 

x P(Na, Br\H E ) dNa dBr 

(2) 

P{Na obs ,Br obs \Na,Br,H L ) 

OO 

x P{Na, Br\H L ) dNa dBr, 

(3) 

where Na obs and Br ons are the Na I and Br 7 equivalent 
width measurements for that star. The likelihood func- 
tions (e.g. P(Na obs , Br obs \Na, Br, He)) are assumed to 
be the product of two independent probability distribu- 
tions, one for Na I and one for Br 7, 

P(Na obs , Br obs \Na, Br, H E ) =P{Na obs \Na, H E ) 

x P(Br obs \Br,H E ). 

(4) 

Each of these terms is modeled as a Gaussian with the 
observed value as the mean and the error in the observed 
value as a obs . For example, the likelihood for measuring 
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Figure 3. Spectral identification of old (red) and young (blue) stars on an LGS AO image of the region. The dashed lines are the outlines 
of each OSIRIS pointing with field name in each corner. The diamond marks the location of Sgr A*. This is the sample of manually 
spectral-typed sources. 



Table 3 

GCOWS spectroscopic completeness 



Mag. Bin a 


E2-l b 


E2-2 


E2-3 


E3-1 


E3-2 


E3-3 


E4-1 


E4-2 


E4-3 


S2-1 


9.5-10.0 






















10.0-10.5 








1.00 












1.00 


10.5-11.0 


1.00 






1.00 






1.00 




1.00 




11.0-11.5 


1.00 


















1.00 


11.5-12.0 






1.00 








1.00 








12.0-12.5 


1.00 


1.00 


1.00 














0.00 


12.5-13.0 


1.00 




1.00 


1.00 




1.00 




1.00 




1.00 


13.0-13.5 


1.00 


1.00 


1.00 




1.00 


1.00 


1.00 


1.00 


1.00 


0.50 


13.5-14.0 


1.00 


1.00 


1.00 


0.50 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


14.0-14.5 


0.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.67 


1.00 


1.00 


0.60 


14.5-15.0 


0.25 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.67 


0.67 


15.0-15.5 


0.08 


0.62 


0.92 


0.62 


1.00 


0.40 


0.67 


0.62 


0.62 


0.00 


15.5-16.0 


0.00 


0.00 


1.00 


1.00 


0.50 


0.33 




1.00 


0.67 


0.00 



a The observed K magnitude bin (not extinction corrected). 

k The fraction of stars that have been manually spectral-typed out of all sources detected in 
NIRC2 imaging in the given magnitude bin. 
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Figure 4. The distributions of equivalent widths for the fainter (K 1 > 14) sample of manually typed early (blue, solid) and late-type 
(red, dashed) stars, along with the best fit Gaussian distribution for Na I (left) and Br 7 (r ight). These distributions arc used as priors for 
the Bayesian inference method of calculating the spectral type probabilities in Section 11.31 
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Figure 5. Completeness as a function of brightness (left) and distance from Sgr A* (right). The spectroscopic completeness based 
on manual typing (solid green) relative to imaging. The imaging completeness based on star planting (dashed blue) are also shown. 
Completeness drops below 50% at K^^ = 15.5; thus, we only analyze luminosity functions and radial profiles down to this magnitude limit. 
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Na EW is: 



P(Na obs \Na,H E ) 



given the measured Bayes Factor, is: 



2lT<TNa 



■ exp 



-(Na-Na obs ) 



(5) 

The likelihood functions in this case are not dependent 
on whether the star is early-type or late-type, since it is 
only a function of our measurements. The priors on the 
distributions of Na and Br, however, are dependent on 
the hypothesis. They are also factored into two indepen- 
dent terms, 



P(Na,Br\H L ) 
P(Na,Br\H E ) 



= P(Na\H L )P(Br\H L ) (6) 
P(Na\H E )P(Br\H E ). (7) 



These priors are modeled as Gaussian distributions based 
on the equivalent width measurements of stars with man- 
ually determined spectral types described in Section [4.1 
(Tabic [21 Figure 0}. For example: 



P(Na\H E ) 



1 



27TCT 



■ exp 



JVa„ 



-(Na - Naprior) 2 



(8) 

If the observed value of the equivalent width is far from 
the peak in the prior, the resulting integral in the evi- 
dence will be small, lending less evidence for this hypoth- 
esis. 

The ratio of the evidences for the two hypotheses is 
the Bayes factor: 



BF 



P(Na obsi Br obs \H L ) 
P(Na obs ,Br obs \H E ) 



(9) 



The use of Bayes factors for evaluating the strength of 
the evidence was first advocated by Jeffcreys (1961), and 
has subsequently been used in co smology to ev aluate dif- 
ferent cosmological models (e.g. iTrottal 1200 8]) . A large 
Bayes factor means that hypothesis H E is preferred over 
hypothesis H E , whereas a small value of BF would mean 
that hypothesis H E is preferred over H E . See Figure[S]for 
a simplified example of evaluating the relative strength 
of the evidence (i.e. the Bayes factor) for a late-type 
star versus an early-type star using only the measured 
Na I equivalent width. Figure [8] shows the distribution 
of In BF 



4.3.2. Deriving Probabilities 

The observed Bayes Factor for each untyped star is 
converted into a probability of being early- type (P E ) or 
late- type (Pl) by running a series of Monte Carlo simula- 
tions. These simulations are used to calibrate the effects 
of local noise sources, such as mis-subtraction of back- 
ground Br 7 gas or halos of nearby bright stars. For 
each untyped star in our sample, we simulate and plant 
100 late-type and 100 early-type stars nearby the source 
as described in Appendix |DJ For every planted star, we 
extract Br 7 and Na I equivalent widths (Appendix [Bj 
to compute the Bayes Factors to calibrate how the local 
environment affect them. 

For each untyped star in our sample, the probabil- 
ity that it is early-type, P E (BF), or late-type, P E (BF), 



Pe(BF) 
Pl(BF) -- 



IeTIe + IlTIl ' 



(10) 
(11) 



where f E is the fraction of simulated young stars with 
Bayes Factors greater than the measured Bayes Factor 
of the untyped source, and f E is the fraction of simu- 
lated old stars with Bayes Factors less than the measured 
Bayes Factor of the untyped source. Figure [7] shows an 
example of how f E and f E are derived, and H E are 
the prior probabilities of observing a late-type or early- 
type star, respectively (see Section H.3.3l for derivation of 
these priors). Statistical uncertainties on the probabil- 
ities are typically < 5%, assuming Poisson errors in f E 
and /l and errors in the priors. 

In Table HI we present the location and properties of 
the untyped sources, as well as their associated Bayes 
Factor and P E and P E values from equations [TU] and QTJ 
Figure [9] shows the probability of being an early-type 
star as a function of In BF for the sources with manual 
spectral types as well as the untyped sources. 

4.3.3. The choice of priors 

A key component to assigning probabilities are the 
prior assumptions on the intrinsic distribution of early 
and late-type sources within the untyped sample. In 
this case, the prior is the relative probability that a star 
is early or late-type. One choice for the priors would 
be to assign equal probability to a star being early-type 
(IIe) or late-type (n^) to each untyped source. Another 
choice would be to assume that the untyped stars have 
the same relative fraction of early and late-type stars as 
those already spectral-typed in a given magnitude bin. 
Neither of these choices is entirely satisfactory; the first 
choice assumes that we have no information about the 
sources that have not been spectral-typed. The second 
assumes that our sensitivity to the two types of stars is 
the same and that there are no location or magnitude de- 
pendent systematics. For example, in general it is more 
difficult to detect Br 7 than Na I at faint magnitudes 
given the smaller Br 7 equivalent widths and complica- 
tions of background subtraction. Additionally, the early 
and late-type sta rs have radial density profiles with very 
different slopes (iBuchholz et al.l [200a IDo et all l2009bl: 
iBartko et al.1 l2010f h To account for this we factor our 
priors into two terms, 

n^; = H Et sens ^E,R (12) 
nL = likens ILr,,R- (13) 

The first contains the differences in line sensitivity to Na 
I and Br 7, and the second incorporates a radial depen- 
dence in the relative number of early-type to late-type 
stars. 

4.3.4. The line sensitivity prior 

To determine the relative sensitivity to the two types 
of stars, we turn to the results of the star planting sim- 
ulations. We must set a threshold in Bayes Factor for 
calling a star early-type versus late-type. In this way, 
we can then examine the simulations to determine what 
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Figure 6. An example of calculating the evidence for a late-type star compared to an early-type star. The solid line in the two plots is 
the measured Na I equivalent width for this star, modeled as a Gaussian with a mean of 3.0 A, and a = 1.0 A(this is the likelihood for 
observing the Na I equivalent width). The priors for both the late- type (dotted red line) and early- type (dotted blue line) stars are also 
shown. The evidence for the hypothesis that the star is late-type or early-type is the integral of the product of the likelihood and the prior. 
The evidence is the area shown in hashed red and blue, respectively. The Baycs Factor is then the relative area of the two distributions. 
For this example, as expected, the evidence is quite strong that the star is late-type, with ln(BF) m 2. 



fraction should have been detected as either type of star. 
Based on the distribution of Bayes Factor for the man- 
ually spectral- typed sources, we chose the Bayes Fac- 
tor threshold to be ln(BF) > —1.76 in order to be de- 
clared a late-type source and \n(BF) < —2.84 to be de- 
clared an early-type source. These thresholds are cho- 
sen such that 90% of the late-type stars have \n(BF) 
above the late-type threshold and 90% of the early-type 
stars are below the early-type threshold in the manually 
spectral- typed sample (Figure [TU]). Using these thresh- 
olds, we examine the entire set of simulated stars and 
compare the BF distributions with the thresholds. If we 
were equally sensitive to early-type and late-type stars 
(HE,sens = ^L,sens = 0.5), the number of planted late- 
type and early-type stars with BF outside the threshold 
would be equal. While the distribution of Baycs Factor 
for the simulated late- type stars are very similar to those 
that have been manually spectral-typed, the distribution 
for simulated early-type stars are skewed closer to zero 
than the observed distribution (Figure ITOj) . The rela- 
tive fractions of undetected sources are 60% early-type 
stars and 40% late- type stars; this means the late- type 
stars are about 1.5 times more likely to be detected than 
early-type stars in the magnitude range of the untyped 
sources. We therefore use as our priors TL E sens = 0.60 
and T\-L,sens = 0.40. We note that this prior is used as 
the starting point for the Bayes Factor analysis - the ac- 
tual evidence for each untyped source also plays a role 
in calculating the posterior probability that a source is 
early or late-type through equations [10] and [TT1 

4.3.5. The radial distance prior 

In addition to the overall sensitivity of our survey to 
the two types of stars, we also consider the different ra- 
dial profiles of the two populations. For example, a ran- 
domly selected star at a projected distance of 10" from 
Sgr A* has a much higher probability of being a late-type 
star compared to one that is at 1". The distance prior 
is the fraction of early to late-type stars as a function of 



radius: 



IL RtE (R)=£(R)/(l + (;(R)) 
n RiL (R) = l/(l + Z(R)), 



(14) 
(15) 

where, £(i?) is the ratio of the radial surface density pro- 
files: 

scrXe a e r- t - 

£.\ R ) = = a n-r, • ( 16 ) 



Z(R)i 



A L R- r - 



~Rr,e(R) and IlR^i?) are determined iteratively us- 
ing the observations. Initially, we use as the priors, 
the surface density profiles of the manually typed sam- 
ple of stars. This initial prior is not dependent on the 
probabilities from the simulations. We derive probabil- 
ities for the untyped stars using this prior along with 
the line sensitivity prior {U E , se ns = 0.6,LI L:Sens = 0.4). 
Radial surface density profiles are then recalculated us- 
ing the complete sample of stars, which we will use as 
our final radial distance prior. See Appendix [E] for de- 
tails of the surface density profile measurements. The 
resulting power-law slope parameters: A E = 2.6 stars 
arcsec -2 , T E = 0.86 ± 0.13, Al = 2.5 stars arcscc -2 , 
= 0.01 ± 0.14, are used as the radial distance prior. 
We find that the resulting surface density profiles are 
relatively insensitive to whether they are measured with 
the iterative approach or stopping after the initial step; 
the difference in the prior probabilities between the two 
steps are less than 3%. 



5. RESULTS 



5.1. 



K' luminosity function 

K' luminosity functions are constructed for both the 
early and late-type stars. First, the manually-typed 
sample is used alone to construct these distributions by 
summing the number of early-type or late-type stars in 
each magnitude bin. We will refer to these as observed 
distributions, since they are equivalent to the observed 
KLFs reported in earlier works before correcting for in- 
completeness. Second, completeness-corrected KLFs are 
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Figure 7. Examples of the results of star planting simulations for stars with a high probability of being late-type (Top, S3-21) or early-type 
(Bottom, S3-348) are shown with their distribution of Bayes factors (we use the In BF) measured from simulated sources planted near these 
stars. The plots show how the distribution of measured Bayes factor for the planted early-type stars (blue, dashed line) are separated from 
those of planted late- type stars (red, solid line). These distributions allow us to calibrate the measured Bayes factor (black, dotted line) 
for each source by sampling how likely it would be to observe a given Bayes factor if the star is early-type compared to being late-type. 
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Figure 8. The distribution of Bayes factors for stars with man- 
ual spectral-types (early-type: dashed blue, late-type: solid red) 
compared to that of untyped population (dotted black). 



constructed by combining the manually-typed and sta- 
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Figure 9. The relationship between measured Bayes factor 
(In BF), and the probability that the star is young, for the set 
of stars without manual spectral-types. This relationship is not 
perfectly correlated because there is environmental variability be- 
tween the location of different stars as well as differences in priors 
on the spectral- types. 

tistically typed samples (e.g. all stars detected in our 
deep images) and applying a small correction for imag- 
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Figure 10. Top: The choice of Bayes Factor threshold for defin- 
ing early-type and late-type stars is based on the cumulative distri- 
bution of Bayes Factor for the manually classified late-type (solid) 
and early- type (dotted) sources. We chose the Bayes Factor thresh- 
old to be ln(BF) > —1.76 (solid vertical) in order to be declared 
a late-type source and ln(iJF) < —2.84 (dotted vertical) to be 
declared an early-type source. These thresholds are chosen such 
that 90% of the late-type stars have ln(BF) above the late-type 
threshold and 90% of the early-type stars are below the early-type 
threshold in the manually spectral-typed sample. Note that 1 - 
CDF of the early-type stars (blue) is plotted to more easily il- 
lustrate how the two populations are separated. Middle: The 
simulated cumulative distributions are also shown for early-type 
(blue) and late-type (red) stars. The relative fractions of simu- 
lated sources that fall outside the thresholds arc 60% for early-type 
stars and 40% for late-type stars. Bottom: The distribution of 
manually-typed early and late-type Bayes factors along with the 
corresponding simulated sources. 

ing incompleteness. In each magnitude bin, the number 
of early-type or late-type stars is given by the sum of the 
probabilities, divided by the imaging completeness: 

N obs 

NE.rnag = 2^ ^Ea/Cmag (17) 

i 

N obs 

NL,mag = ^ PL,i/C ma g (18) 

i 

where N }, s is the total number of stars, C mag is the 
imaging completeness within the given magnitude bin, 
and Pe,% an d Pl,i are the the probability that the star is 
early or late-type, respectively. The manually spectral- 
typed early and late-type stars are respectively assigned 



either Pe,% = 1 or Pl j = 1. Note that the imaging com- 
pleteness correction is identical for both early-type and 
late-type samples. NE,mag and Nr,^ mag are assumed to 
have Poisson errors^. Poisson errors on the final num- 
ber are a good approximation of the true error since we 
know the number of stars very well from deep images 
(i.e., imaging completeness is very high), and the statis- 
tical uncertainties in the probabilities for early and late- 
type classification are assumed to be negligible. Figure 
[TT1 shows both the observed and completeness-corrected 
K' luminosity function for the early and late-type stars. 
In the faintest bin (15.0 — 15.5 mag), 96% of the stars are 
detected in imaging, of which about 50% of those have 
manual spectral- types. The early-type KLF increases 
smoothly with fainter magnitudes, while the late-type 
stars exhibit a large jump in the faintest bin due to the 
presence of red clump stars. 

5.1.1. Luminosity function of the S-stars 

We also compare the luminosity function of the S-stars, 
defined here as early-type stars wit h R < 1.0" (similar to 
the definition in lGenzel et al.ll2010T > . to that of early-type 
stars further out. We split the sample of early-type into 
different regions: (1) R < 1", (2) 1" < R < 12", and 
(3) all R. Figure Q2] shows the completeness-corrected 
luminosity functions of these three samples of stars. The 
stars with R < 1" arc 100% complete to an extinction 
corrected K' AA < 15.5 (11 stars: SO-2, S0-1, SO-3, SO-5, 
S0-11, S0-4, SO-9, SO-31, S0-14, Sl-3, SO-15). This high 
completeness is due to the lower than average extinction 
in this region as well as deeper spectro scopic observa- 
tions . As observed by other studies (e.g. IPaumard et al.l 
|2006[) . the central 1" has a lower density of stars with 
K' AA < 14.0, compared to the outer region. However, 
for the stars with K' AA > 14.0 in the central arcsecond, 
the luminosity function is statistically consistent with the 
faint end of the luminosity function for stars at R > 1", 
and for the total population. 

5.2. Radial surface- density profiles 

We use a Bayesian parameter estimation method to 
determine the surface-density profiles of the early and 
late-type stars. Previous observations of the surface- 
density profiles of stars in this region were estimated by 
fitting a power-law to number counts binned by radius. 
The process of binning can lead to large variances in 
the inferred power-law fits, especially if the population 
is to be separated by luminosity, for example. In order 
to overcome some of these limitations, we use an un- 
binned fit to a power-law model, £ cx R~ T . We compute 
the Bayesian posterior probability distribution for the 
power-law slope T, using the individual star's positions 
as the likelihood, and assuming a flat prior for T. See Ap- 
pendix [E] for more details. The power-law slopes of the 
observed (as defined in Section 15. ip surface-density pro- 
files for the early and late-type stars is Te — 0.93 ± 0.11 
and Tl = — 0.04± 0.10, respectively. We plot these pro- 
files in Figure 1131 along with radially binned points for 
illustration and comparison. The completeness-corrected 
late-type surface-density profile has a best fit power-law 

10 In cases where N = (e.g. crjv = vN), we conservatively 
adopt opf = 1. 
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Figure 11. The luminosity function of the early (blue, solid) and late-type (red, dashed) stars binned in 0.5 mag bins. Left: the observed 
luminosity function using the sample of manually typed stars and the observed K' magnitudes. Right: the differential-extinction and 
completeness-corrected luminosity function. 

sequence (MS) stars and the more massive O stars and 
OB supergiants. This split is also motivated by the 
fact that there appears to be a significant lack of young 
stars brighter than this threshold in the central 0."8, 
which has been noted by many previous observers (e.g. 
iPaumard et all 120061 : iBartko et all 12010( 1 . We wish to 
investigate whether there is a difference in the density 
structure between the B MS populations and the more 
massive young stars. We find that while the fainter B 
stars have a marginally shallower density profile, with 
r faint = 1.06 ± 0.25, than the brighter population, with 
r bri ght = 1.25 ± 0.22 (for R > 1.0"), their profiles are 
statistically consistent with a single power-law for pro- 
jected distance R > 1" from Sgr A*. When including 
the area inside of 1", the power-law fit for all faint B- 
type stars has T = 1.06 ± 0.13, consistent with the fit 
for R > 1". Figure fl4l shows the profiles of the different 
samples of stars. It is also unclear at this point whether 
a single power-law is a good fit to the density profiles 
of both populations, as there appears to be a plateau in 
the surface-density profile from about 1" to 4", beyond 
which the surface-density drops. This may indicate that 
a broken power-law is a better fit to the surface-density 
profile. However, because of the small number of stars, 
this effect is not statistically significant at this time. 
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Figure 12. The completeness-corrected luminosity functions of 
the early-type stars as a function of projected distance from Sgr 
A* are shown. Blue, solid - all early-type stars in the survey. Or- 
ange, dotted - early-type with projected distance of 1" < R < 12" 
from Sgr A*. Black, dashed - early-type stars within 1.0" of Sgr 
A*, scaled by a factor of 0.18 to better compare the slope of the lu- 
minosity function. The bright end (K^ A < 14.0) of the luminosity 
function for the central 1.0" is inconsistent with that found further 
out, lacking bright stars. However, the faint end (K'^ A > 14.0, 
B-stars) of the S-stars luminosity function is consistent with the 
luminosity function for stars at > 1". 

of = 0.16 ± 0.07, while the early-type surface-density 
profile has a best fit slope of T E = 0.93 ± 0.09. For 
comparison, a traditional least-squares fit to a binned 
radial profile, with each radial bin containing roughly 
equal number of stars, has Tl = 0.12 ± 0.16, and 
Te = 0.83 ±0.14. While the two methods are consistent, 
the Bayesian method has the advantage of utilizing the 
precise positions of stars rather than binned positions, 
resulting in smaller uncertainties. Table |4] summarizes 
the fits to the different populations of stars. 

We also examine the early-type population for evidence 
of mass segregation by examining the surface-density 
profiles for stars brighter than K' AA = 14.3, compared 
to the fainter population. This cut is chosen because it 
is approximately at the division between B-type main 



6. DISCUSSION 

6.1. K' luminosity function: early-type young stars 

We find our measured luminosity function for young 
stars between 0."8 < R < 12" fr om Sgr A* to be mu ch 
steeper than the one reported by IBartko et alJ (I2010D in 
the sa me region. The luminosity function in lBartko et al.l 
(|2010( 1 is essentially flat between K s = 12 and K s = 16. 
In comparison, the luminosity function in this study rises 
continuously toward fainter magnitudes. Figure [T5l com- 
pares the extinction and c ompleteness-corrected luminos- 
ity of lBartko e t al.l (|2010f ) to the present study for early- 
type stars at a projected distance of 0."8 < R < 12" 
binned in 1 magnitude bins. The luminosity functions 
are scaled to have the same value at K' = 12.0 in or- 
der to compare the differences in the slopes. Even be- 
fore completeness correction, our survey has a compara- 
ble fraction of young stars in the 14.5 < K' < 15.5 to 
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Figure 13. The azimuthally averaged surface-density profile of early-type (blue) and late-type stars (red). The most probable surface- 
density power-law slopes are also plotted (early: dashed, late: solid). Left: The number counts have not been corrected for completeness. 
These consist of Sample 1 stars. The most probable fit for the two populations are: Tl = —0.04 ± 0.10, and Te = 0.93 ± 0.11. Right: 
completeness-corrected surface-density profile (Sample 4). The most p roba ble fit for the two populations are: Fl = 0.16 ± 0.07, and 
r^; = 0.93 ± 0.09. For details of the completeness correction, see Section 14.31 

fcrred mass function compared to lBartko et al.l (|2010D . It 

is not trivial to derive an IMF from a luminosity function 
as there are many variables that can affect the luminos- 
ity besides mass. For example, the age, star formation 
history, and metallicity will all affect the tra nsformation 
from mass to luminosity. In ILu et all (|2012|) , a detailed 
analysis is performed using the data presented in this pa- 
per and employing a combination of stellar evolution and 
stellar atmosphere models. Here, we will focus our dis- 
cussion on a direct comparison of the observed luminosity 
functions and our approach to completeness correction. 

One possible explanation for the discrepancy between 
our KLF and previous measurements is our different 
method of correcting for spectroscopic completeness. 
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Figure 14. The azimuthally averaged, completeness-corrected 
surface-density profile of young stars (blue). Also, the surface- 
density profile of those that are brighter (orange) and fainter (pur- 
ple) than K' = 14.3, approximately the magnitude that separates 
B MS stars from the brighter OB supcrgiants and WR stars. There 
is a large drop in the density of bright stars in the center, but the 
outer radial profiles are consistent between the two populations. 
Power law slopes and errors for the various populations are given 
in Table IH 

that of the completenes s-corrected luminosity function 
from lBartko et all (|2010f ). Any completeness corrections 
would then increase the steepness of the K 1 luminosity 
function in this study. The steeper slope of the luminos- 
ity function found in this study will result in a steeper in- 



Our method takes advantage of our nearly complete 
knowledge of the location and brightness of stars from 
deep imaging; only the spectral-types of some stars are 
unknown. We incorporate this knowledge, along with in- 
formation in the spectra of these untyped stars, to assign 
a statistical probability of being early-type or late-type. 
In comparison, traditional completeness corrections ig- 
nore all information on un-typed stars and only plant 
simulated stars to estimate a completeness correction to 
be appli ed to the manually- typed population. Further- 
more, in lBartko et al.l (|2010t) . the simulated stars appear 
to be typed in a manner that is different from the ob- 
served stars; the simulated stars are declared early-type 
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Figure 15. A comparison between the K'^ A luminosity function 
reported in this paper and that reported in the literature. The 
plot compares the observed (purple asterisks) and completeness- 
corrected K'^ A luminosity function of early-type stars (blue solid) 
located at a projected distance 0."8 < R < 12" from Sgr 
A*, with the comp leteness- corrected K luminosity function from 
IBartko et al.l Q2010T ) (black diamond) in the same radial range. 
The numb er counts are normal ized between our observations and 
those from IBartko et alj H20101 ) in the K' AA = 11.5 - 12.5 bin in 
order to compare the relative difference in the slopes (the error 
bars are also scaled accordingly). The bin centers for the ob- 
served K Ayl luminosity function are shifted slightly for clarity. We 
find a signifi c antly greater fraction of faint young stars than in 
IBartko et al.l 1 120101) . The observed luminosity function is 

near ly identical to the co mpleteness-corrected luminosity function 
from|B"artko ct al. (20fi]), so any amount of incompleteness correc- 
tion will lead to a greater number of faint B stars. 



based on CO alone, while the observed stars are de- 
clared early-type when they have Br-gamma or He I, and 
lack CO lines. This may lead to an over-estimate of the 
completeness to early-type stars, in which case, their re- 
ported numbers of early-type stars are underestimated in 
the faint bins where the completeness corrections domi- 
nate. 

Anot her significant differ ence between this survey and 
that of IBartko et alj (|2010() is the region covered by the 
two surveys; this survey is done predominantly in the 
direction of the proj ected clockwi s e you ng stellar disk, 
while the survey from lBartko et al.l (pOlOt) covers a region 
largely perpendicular to the disk plane (Figure [2]). The 
present survey is likely to contain more star s that belong 
to the clockwise disk of young stars, while IBartko et alj 
(|2010f ) showed that only one star fainter than K = 15.0 is 
consistent with being on this disk. A detailed kinematic 
analysis of the current survey is necessary to place similar 
constraints on disk membership (Yelda et al. in prep.), 
but the differences in the luminosity functions, if it is not 
due to differences in completeness corrections, may indi- 
cate a difference in IMF between those stars on the disk 
and the field population represented by the observed B 
stars. This would be the first indication that the stellar 
population is different between these two populations; 
previous observations of the OB supcrgiants and O main 
sequence stars have shown indistinguishable differences 
in number or age of those stars. A difference in the num- 
bers of lower-mass main sequence stars could indicate 
differences in their formation or in their dynamical evo- 
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Figure 16. Comparison between the completeness-corrected and 
dereddened luminosity function of the late-type giants observed 
with OSIRIS (solid, black) and the K lum i nosity function ob- 
served in Baade's Window from Ticdc et al. (1995), which gives 
a deep sampling of a bulge field near the Galactic center. The 
OSIRIS luminosity function is cut off beyond K = 13.0 (observed 
K' ~ 15.5), whe re our completenes s correction is less reliable. The 
amplitude of thejTiede ct al. ( 1995) luminosity functions are scaled 
to the OSIRIS luminosity function to match at the K = 11.0—11.5 
magnitude bin. 

lution. With the currently published data sets, it is not 
yet possible to quantitatively assess either of these sce- 
narios. A larger systematic survey of the B stars, along 
with their kinematics to establish disk memberships, will 
be necessary to address these questions. 

6.2. K' luminosity function: late-type giants 

We find that the late-type luminosity function in this 
study is very similar to that of the inner bulge popula- 
tion. The late-type luminosity function is comparable to 
that of Baade's Window, a field located about 4° from 
the Galactic center with very low extinction (Ak ~ 0.14) 
and so has been studied extensively in the past. Fig- 
ure 1 161 shows the K luminosity function from lTiede et all 
(|1995l ). who combined dee p observations of Baade's Win- 
dow with ear l ier wo rk by iFrogel fc Whitfordl (|1987l ) and 
iDepoy et al.1 |l993). These observations are corrected 
for reddening by dust, but not by the distance modu- 
lus (K Q ), which reaches a depth of K Q = 16.5, corre- 
sponding to about K — 19.2 at the Galactic center when 
the additional extinction is included. We deredden the 
K' luminosity function from the OSIRIS spectroscopic 
sample in order to compare with the one from Baade's 
Window. We scale the amplitude of the luminosity func- 
tion of Baade's Window to the K Q = 11.0 — 11.5 bin 
of the completeness and dereddened luminosity function. 
We find that the two generally agree well, down to the 
spectroscopic limit at K' = 15.5, or K a = 12.5 — 13.15, 
where the red clump stars are concentrated in the lumi- 
nosity function. Below K' = 15.5, we no longer have 
spectroscopic differentiation between the young and old 
population, but it is clear from Figure [T5] that the GC 
luminosity function matches the red clump features from 
Baade's Window. The match in the slope and location of 
the red clump to that of the bulge indicates that the star 
formatio n history at the Ga l actic center may be similar 
(see also lManess"et^ [20071: TPfuhl et al.|[201ll) . 
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6.3. Surface- density profile of young stars 

The power-law fit to the projected surface-density pro- 
file of all the young stars is consistent w ith that of pre- 
vious observations (T = 0.90 ± 0.09). iPaumard etaTl 
(|2006f ) reported a surface-density profile power-law slope 
of r = 2.1 ± 0.17 in the plane of the young stellar disk, 
including only stars in the clockwise disk. As the current 
study does not separate the young stars into disk mem- 
bers and non members, we must compare the surface- 
density profiles for all y oung stars. Us i ng th e entire 
sample of young stars in IPaumard et al.1 (|2006|) within 
the field of view of this study, regardless of disk mem- 
bership, the best fit power-law to the projected surface- 
density profile is T = 1.12 ± 0.13, consistent with the 
measureme nt here. Our results ar e also consistent with 
those from IBuchholz et al.1 (120091 ) . who used medium- 
band imaging to separate early-type stars with low CO 
equivalent widths from la te-type stars wi t h stro ng CO 
absorption at ~ 2.3 /xm; IBuchholz et~aTl (|2009j ) found 
that at a projected distance of 1" to 10" from Sgr A*, 
the early-type stars have V = 1.08±0.12. With the iden- 
tification of these young stars, further insights into their 
origins can be obtained with the addition of kinematic 
data to reconstruct their orbital properties (Yelda et al. 
in prep.). 

6.4. Properties of the S-stars 

Our observations show that other than a deficit of stars 
brighter than K' AA < 14.0, the early-type S-stars within 
the central ~ 1" have very similar properties to those 
found at greater distances. There are no significant dif- 
ferences in the surface-density profile of the B-type stars 
(stars with < 14.3) within and outside the cen- 

tral arcsecond (Table [4]). The B-stars show a contin- 
uous surface density profile throughout the survey re- 
gion. The luminosity function within the central arc- 
second is also consistent with the early-type stars fur- 
ther out for K' A A > 14.0 (Figure [TSJ) - In contrast, 
IBartko et al.l (|2010| ) found that the luminosity function 
of stars with R < 0."8 is significantly steeper than that of 
stars with 0."8 < R < 12". For comparison, we also con- 
structed a luminosity function for early-type stars with 
R < 0." 8. The lumin o sity fu nction in this region is the 
same as IBartko et al.l (|2010f ), because of the high com- 
pleteness in this region for both surveys (the spectral 
ide ntifications are the same in this region). We find, as 
in (|Bartko et al.l[20Tot ). that this region is missing bright 
stars with K' AA < 14.0 compared to the luminosity func- 
tion at R > 0.8". However, we find that the faint end 
of luminosity function (K' AA > 14.0) inside R < 0."8 
is consistent with the faint end of the luminosity func- 
tion for stars with 0."8 < R < 12", due to the steeper 
measured luminosity function in this survey. Our results 
suggest that the B-type S-stars may represent a continu- 
ous population of B-stars throughout the central 0.5 pc. 
It is unclear at this time however, whether the S-stars 
can have originated from the most recent star formation 
event that formed the young disk of stars further out; for 
example, an explanation for the curious deficit of bright 
stars in the inner 1" will be necessary for this hypothe- 
sis. Resolving this issue will have a strong impact on our 
understanding of the star formation in the region as well 
as the timescales for the dynamical mechanisms that are 



necessary to bring the S-stars so close to the supcrmas- 
sive black hole. 

6.5. Cusp clearing out to 0.5 pc 

The observed flat surface-density profile of the old red 
giants extends out to the edge of our survey at about 0.5 
pc, about a fac tor of 3 further t han our initial spectro- 
scopic survey in lDo et al . (2009a). While this survey pre- 
dominantly samples the region east of Sgr A*, it should 
be representative of the distribution of old stars in this 
region as there do not appear to be any detectable devia- 
tions from spherical symmetry. This large core profile is 
also consistent with th e narrowband imaging results from 
Buchholz et al.l (120091) and the spectroscopic results from 
Bartko et all (|2010f ) based on samples more to the north 
of Sgr A*. It is unclear from the present data whether 
there is a break in the surface-density profile at larger 
radii, as the survey truncates at about 12". Using a bro- 
ken power-law model, we have attempted to constrain 
the break radius and the outer power-law slope, but the 
current data have insufficient radial coverage to strongly 
constrain the location of the break. Because the pro- 
jected surface density profile is so flat, it is also difficult 
to determine the true spatial density profile of the late- 
type stars (p(r) cx r -7 , where r is a 3D distance). This 
limitation can be removed with the inclusion of kinematic 
information. For example, Jeans modeling of proper mo- 
tion and rad ial velocity mea surements has successfully 
been used bv lDo et" al. (2012]) to constrain the power-law 
exponent, 7, within the survey region presented here. 

7. CONCLUSIONS 

In this paper, we report the results of our new spec- 
troscopic survey of the central 0.5 pc of the Milky Way. 
This study presents both new observations as well as new 
methodologies. Our new data extend previous spectro- 
scopic observations along the disk plane of the young 
stars by a factor of ~ 2 — 3. We develop a new method for 
statistical spectral-typing that we use for completeness 
correction. This method allows us to takes advantage 
of prior information about the stars. Most importantly, 
this includes information about the locations, bright- 
nesses, and spectra. Because over 95% of stars within 
the magnitude range of interest are detected in imag- 
ing, this provides a very robust method for constructing 
completeness-corrected luminosity functions and surface- 
density profiles for young, early-type stars and old, late- 
type giants. 

We find that the measured radial surf ace-density pro- 
files are consistent with previous studies fB uchhol zet al.1 
12001 IDo et al.ll2009al ; IBartko et al.ll2010D . The surface 
density profile of the late-type stars appears flat within 
our survey region, suggesting that the 'core' in the red 
giants is at least ~ 0.5 pc in size. The early-type stars 
have a much steeper radial surface density profile such 
that they dominate the stellar density within < 0.04 pc 
from Sgr A*. 

The luminosity functions of both the late and early- 
type stars rise towards fainter magnitude bins. The late- 
type stellar luminosity function is consistent with the 
inner bulge of the Galaxy, indicating that the star for- 
mation history of the Galactic center may be similar. 
The luminosity function for the early-type stars is con- 
sistent with that of previous studies at the bright end 
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(K' < 13.0; iPaumard et al.|[200l iBuchholz et al.ll2009D . 

but is ste eper than reported in a recent study at the 
faint end (|Bartko et al.l I2010T) . This steepening of the 
faint end of the luminosity function will likely result in a 
steeper mass function t han t he very top-heavy IMF re- 
ported in IBartko et al.l (|2010|) . The derivation of a mass 
fu nction from a lu minosity function is presented in Paper 
II (ILu et al.H20Ll . 

We find that the S-stars at R < 1" and K' AA > 14.0 
(B-type) have the same luminosity function and surface 
density profile as the B-type stars further out. This sug- 
gests that the population of all B-type stars in the central 
0.5 pc may be related, though there is insufficient infor- 
mation at this time to determine whether they originate 
from the same star for mation event as the young stellar 
disk (see also, Paper II. ILu et al.ll2012h . 

Accurate measurements of the luminosity function are 
important, as different luminosity functions can lead to 
very different conclusions about star formation in the ex- 
treme tidal environment of the Galactic center, which af- 
fects our understanding of star formation in general. It is 
important to note that current spectroscopic studies are 
limited to K' < 15.5 mag , which corresponds to ~ 10M@ 
(Paper II ILu et aT1l2012[ ). To more completely compare 
the Galactic center mass function with local star-forming 
regions, observations down to about a solar mass are nec- 
essary. This mass corresponds to about K' sa 21 mag at 
the Galactic center, which cannot be reached by cur- 
rent IFU instruments. This regime for scientific study 
will only be opened with future Giant Segmented Mirror 
Telescopes (GSMT), such as the Thirty Meter Telescope 
(TMT). 
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Table 5 

Stars observed in K broad-band for spectral-type verification 



Name 


K' 


ln(BF) 


Kn3 Sp. Type a 


Kbb Sp. Type 


Field 


S2-317 


15.52 


-2.73 


Late 


Late 


S 


S2-55 


15.21 


-0.68 


Unknown 


Late 


SE 


S2-61 


15.36 


1.00 


Unknown 


Late 


E 


S2-64 


15.57 


-1.13 


Unknown 


Late 


E 


S2-77 


13.38 


0.11 


Late 


Late 


S 


S3-7 


13.56 


-3.50 


Late 


Late 


SE 


S3-96 


14.31 


-1.96 


Unknown 


Early 


sw 


S4-46 


14.72 


1.84 


Unknown 


Late 


E2-2 


S5-127 


15.62 


-0.41 


Unknown 


Late 


E2-3 


S5-211 


13.21 


-1.00 


Unknown 


Late 


E2-3 


S5-237 


13.21 


-2.92 


Unknown 


Early 


E2-1 


S6-77 


14.11 


-3.13 


Late 


Late 


E2-3 



a Initial spcctral-typc using only Kn3 spectra with the criteria described 
in this paper in Section ICT 



APPENDIX 
SPECTRAL-TYPE VERIFICATION 

The identification of spectral types is largely based on the presence or absence of the Br 7 an d Na I lin e s. The are, 
however, a few bright sources for which there are no spectral features or very weak Na I lines. In lDo et al.l (|2009aD . we 
classified the bright, K' < 14 featureless stars within the Kn3 wavelength range as young, as they were p resumed to be 
OB st ars with very low or zero equivalent width in Br 7, consistent with the variance in Br 7 reported bv lHanson et al.1 
(|1996D . However, at the faint end of this range (K' = 13-14), we have determined that a few sources contain very weak 
Na I features, with equivalent widths < 2 A, compared to the average Na I equivalent width of ~ 3.8 A. In order to 
determine whether these sources, and other similar ambiguous sources, are young or old, we obtained K-broad band 
spectra with OSIRIS for a sample of stars having Na I equivalent width < 2 A(Tablc [5]). Figure [T71 shows examples 
of the low Na equivalent width sources compared to more typical early-type sources, and the corresponding spectra 
in the K broad-band filters. The K broad-band spectra cover the CO band-heads at 2.3 ^m, which are very strong in 
late-type stars and are much better discriminators of the temperature of the star. We find that the sources with low, 
but detectable, Na I equivalent widths in the Kn3 filter have detectable CO features, indicating that they are late-type 
stars. However, the CO equivalent widths are smaller than the majority of the K and M giants. This is consistent 
with those of warmer giants of ~ 5000 K instead of with temperatures in the range 3000 - 4000 K (jForster Schreiberl 
I2000D for K and M stars. These stars likely represent a younger population (100-300 Myr old) compared to the 
~ 1 Gyr old M and K giant s, but we will classify them as late- type to separate them from the much younger ~ 6 
Myr old population (see also lBlum et alj|2003t iPfuhl et al.l 120 111) . Based on these observations, we make two minor 
modifications for the spectral- type classification criteria: (1) the brightness cut off for featureless stars to be classified 
as early-type is moved from K' < 14.0 to K' < 13.0, (2) sources with small, detectable Na I equivalent widt hs are 
classified as late-type. This results in the reassignment of the spectral types of 8 stars from IDo et al . (2009a) from 
early-type to late-type, 4 of which have K' > 15.0. We find that these modifications have resulted in a much more 
robust method for spectral classification, with no incorrect assignments when comparing the followup sample between 
Kn3 and Kbb. We include the properties of stars classified using Kbb spectra in the Bayesian inference model through 
the distribution of equivalent widths in Br 7 and Na I, as measured using the Kn3 spectra (Section 14. 3j) . 

MEASUREMENT OF SPECTRAL FEATURES 

In order to facilitate source extraction and identification, and to run Monte Carlo simulations for completeness 
corrections, we also develop a method for automating the measurement of the equivalent widths of the Br 7 and Na I 
doublet lines as well as the radial velocity of the sources. The steps in this procedure are as follows: 

1. We cross-c orrelate the spectru m with that of M3II giant HD40239 from the SPEX telescope infrared spectral 
templates (|Ravner et a.1.1 [200~9h in the range of the observed filter from 2.121 to 2.220 ^m; this template was 
chosen because it has a high peak correlation value when cross-correlated with most of the observed late-type 
sources. 

2. If the correlation coefficient is greater that 0.5, then the radial velocity is measured based on the peak of the 
cross-correlation function. We determine the thresholds for correlation from the K' > 14.0 subset of manually 
spectral- typed sources (Sample 2). The location of the cross-correlation peak is determined by fitting a parabola 
to the five closest points around the pixel with the maximum correlation. Once the radial velocity is determined, 
we then shift the spectrum to rest wavelengths and measure the equivalent width of the Na I do ublet. The doublet 
is mea sured by integrating the continuum-removed spectrum between 2.2053 to 2.2101 fim as in lForster~ Schreiber 

poooh . 
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Figure 17. Left: Spectra of sample sources in the Kn3 filter illustrating the types of stars seen in our sample. From top to bottom, the 
spectra are: 1) a typical late-type stellar spectrum with prominent Na I lines at 2.206 and 2.2090 fim; 2) and 3) are late-type stars with 
small Na I equivalent widths; 4) an early-type star with a featureless spectrum in Kn3; and 5) a typical early-type star with Br 7 absorption 
at 2.1661 (mi. Right: The K broad-band spectra for the same sources, showing how the CO band-head can be used as discriminator of 
early-type versus late-type for the more ambiguous sources in Kn3. Featureless early-type sources in Kn3, such as S3-10, may also have He 
absorption in the broad-band filter. 

3. If the cross-correlation coefficient is below 0.5, the wavelength range over which we performed the cross correlation 
around the Na I doublet is reduced to between 2.20 to 2.215 /im. By restricting the wavelength range, the 
sensitivity to the Na I feature is increased for low SNR spectra. If the peak of the cross-correlation function is 
now greater than 0.5, the radial velocity and the Na I equivalent width is measured in the same way as in step 
2. If the correlation peak is still less than 0.5, we do not apply any velocity shifts to the spectrum, but still 
integrate over the region around Na I doublet to establish the equivalent width within this region. 

4. We then measure the equivalent width around Br 7 by first cross-correlating the spectrum at ±4000 km/s 
around the Br 7 line with a template spectrum constructed from multiple observations of SO-2 and shifted to 
rest wavelengths. If the cross-correlation coefficient is greater than 0.3, we fit a Gaussian to the region ±0.1 /j,m 
around the wavelength corresponding to the peak lag in the cross-correlation function. If the peak correlation is 
less than 0.3, we fit a Gaussian to the wavelength region 2.1661 ±0.1 /im. Because we are interested in measuring 
the faint B stars, which should have their Br 7 in absorption, we restrict the Gaussian fit to absorption features. 
We also require that the width of the Gaussian fit be greater than 1 spectral pixel in order to avoid fitting cosmic 
rays or bad detector pixels. 

5. The errors on the measured equivalent widths and radial velocities are estimated by splitting the data into three 
subsets. In order to obtain comparable SNR between these three spectra, we sort all of the spectra for each 
star by their SNR and populate each subset such that the resulting combined spectra would have similar SNR. 
We then apply steps 1-4 to each of the three spectra. We use the standard deviations of the radial velocity and 
equivalent widths as the error for the corresponding measurement. Most stars have between 6 and 9 spectra 
observed, though there are stars with as few as one measurement if it is at the edge of our dither pattern. For 
stars with fewer than three measurements, we estimate the error using a fit for the correlation between SNR and 
equivalent width uncertainties. Using the power-law fit to this relationship and the given SNR of the spectrum, 
we infer the error on the measured parameters. 

6. To filter out spurious fits, we do not consider stars with Na I or Br 7 equiv alent widths > 15 A. This threshold is 
set a bove what is physically expected for any star at the Galactic center ([Hanson et al.l ll996: Fo rster Schreiberl 
I2000T ). Stars with spurious detections arc flagged and those measurements arc not considered in the subsequent 
analyses. Approximately 13% of the stars in our sample of stars with no manual spectral types (Sample 3) 
have either equivalent width measurements above this threshold using the automated routine, or the routine was 
unable to provide an equivalent width measurement. These stars tend to have have SNR < 5. 

Because these measurements are fully automated, they may be susceptible to systematic errors in the spectra that 
would lead to poor estimates of the equivalent widths of the spectral lines. To test the accuracy of this automated 
method, we compare the results to those of measurements that have been individually extracted and checked by 
eye. For the late-type giants, we find that the measurements using the automated routine are consistent with manual 
measurements; the mean of the distribution of Na I equivalent width for old stars with brightness K' < 14.0 is 3.8 ±1.3 
A compared to 3.9 ±1.2 A measured manually. For the early- type stars, we find the mean equivalent width of Br 7 
for young stars with K' < 14.0 from the automated routines (3.2 ± 1.6 A ) to be consistent with the equivalent widths 
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measured manually (3.9 ± 2.4 A). Sec Figure IT51 for plots of the distributions of equivalent widths of Br 7 and Na I. 
We conclude that the automated and manual procedures achieve comparable measurements. 




Na I Eq. Width (Angstrom) Br y Eq. Width (Angstrom) 

Figure 18. Left: the observed distribution of Na I equivalent widths (solid black) as measured manually compared to that using the 
automated procedure described in Section [4] for sources with K' > 14.0. Right: a similar plot for the distribution of equivalent widths for 
Br 7 at 2.1661 fim. These distributions show that the automated algorithm produces similar results to measurements that require human 
interaction. 



PHOTOMETRIC CALIBRATION 

Star lists are photometrically calibrated using published magnitudes reported in lSchodel et al.l (|2010l ). They provide 
an extensive star list with absolute H, Ks, and L' photometry and we select calibrator stars as those with brightnesses 
of H<18, Ks<16, L'<15, and photometric errors of <jh < 0.03, <tk s < 0.02, and 07/ < 0.04. We also exclude stars 
that have neighbors within 0'.'3 and AKs < —1.5. The Ks magnitudes are then converted into our K' filter set using 
the equation 

K' = Ks + 0.00683 + 0.01049 * (H - Ks) (CI) 

This K'— Ks conversion equation is determined by simulating synthetic spectra for a stellar po pulation with a n age 
of 5 Gyr at a distance of 8 kpc using models of stellar evolution and atmosp heres described in dLu et aT1l2012[ ). The 
synthetic spectra are reddened using the Galactic center extinction law by iNishivama et al.l (|2009f> and a range of 
extinction values (Ak s = 2.4-3.0, AAk s = 0.1). The reddened synthetic spectra are convolved with atmospheric and 
filter transmission profiles using the package pysynphot to generate synthetic photometry for a suite of near-infrared 
filters, including H, K', and Ks. A linear relation is then derived for the subset of simulated stars that fall along the 
red giant branch from Ks=14-17, since most of the observed calibrators are cool, red giants (Figure [19] Simulations of 
a younger (6 Myr) population shows that the K'— Ks conversion for hot, young stars would result in K' photometric 
differences of less th an 0.01 mag. This is far less than the 0.06 mag zeropoint error in the calibrator Ks magnitudes 
(jSchodel et alJ l2010h . The final photometric errors for the sample of stars that are brighter K'=15.5 has a mean of 
0.08 mag. 



IMAGING COMPLETENESS 

The NIRC2 imaging completeness as a function of position and brightness is estimated by planting simulated stars 
and determining how well they can be recovered. To accurately estimate the completeness, the images containing 
simulated stars must be analyzed in exactly the same fashion as the real images. We therefore plant stars in each tile 
of the mosaic as well as each tile's three subset images at the same position and brightness. The simulated images for 
each tile are then analyzed with the same requirement that sources must be detected in both the combined tile image 
and all three of its subset-images to be identified as a star. It is not necessary to plant stars at the same position in 
multiples tiles where the tiles overlap as the mosaic process did not impose additional requirements for sources to be 
detected in more than one tile. 

Artificial stars are generated in a grid in both magnitude and position. Simulated magnitudes range from K'=7.7 — 
19.5 in steps of AK' = 0.25 mag. Simulated positions are set in a regular grid separated by 0'.'25. This grid of artificial 
stars cannot be planted in a single simulated image without dramatically impacting the stellar density and resulting 
completeness measurements. Therefore many si mulated images are g enerated, each one containing artificial stars at a 
fixed brightness and spaced 0'.'5 apart (see also. iSchodel et afll2007D . For a given brightness, 4 simulated images are 
required to achieve the final 0'.'25 spatial sampling. Over 300 simulated images are necessary to cover the full brightness 
range just for the combined image of each tile. In total, 16,784 simulated images were generated for the 13 tiles in 
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Kp-Ks = 0.0068 + 0.0105 * H-Ks 
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Figure 19. The K' - Ks correction as a function of H - Ks color derived from synthetic photometry of red giant stars observable at the 
Galactic center. Synthetic stars are selected from a simulated 5 Gyr starburst with solar metallicity at a distance of 8 kpc and reddened 
using the Nishivama ct al. (2009) Galactic center extinction law. Stars on the red-giant branch and in the red-clump are the most numerous 
in Galactic center observations, which corresponds to selecting stars from the synthetic isochrone that have 14 < Ks < 17 at Ajf s = 2.7. 
These stars are shown for extinctions ranging from 2.4 < A^ s < 3.0, the observed range of extinctions in th e central parsec of the Galaxy. 
A linear relation between K' - Ks and H - Ks is fit and used to convert between Ks magnitudes reported in[Schodel et al. (201(3) and the 
K' magnitudes reported in this paper. 

2006, including their total and three subset images, and were analyzed using Starfindcr in an identical manner to the 
observed images. This required ^1400 computing hours, or about 1 week when run in parallel on several powerful 
desktop computers. The simulated and recovered starlists for all the tiles were mosaicked together in the same manner 
as the observed starlists. The final outcome of the star planting simulations is a 3D cube of the number of simulated 
and recovered stars at different X, Y, and K' values, allowing completeness curves to be calculated for different areas of 
the image. Star planting simulations are extremely time consuming; therefore completeness maps are constructed only 
for a single epoch, 2006, which has slightly higher Strehls than the other two epochs. This may tend to overestimate 
the completeness; however, the effect is negligible at K'<16 where luminosity functions are analyzed in this paper. In 
this work, the completeness is calculated for the entire OSIRIS field of view. Figure [5] shows the resulting completeness 
curve and how completeness changes with radius. Completeness decreases inside of 4" , primarily due to the limited 
contrast around the bright Wolf-Rayet stars concentrated in this central region around Sgr A*. 

OSIRIS STAR PLANTING SIMULATIONS 

Our effort to spectral-type the stars using the OSIRIS Kn3 data are subject to several limitations, which results in 
incompleteness in our survey for stars with K' > 12.5 mag. The following is a list of contributors to the incompleteness 
of the survey: 

• Intrinsic variations in the equivalent widths of absorption lines will make some stars more difficult to detect than 
others at a given SNR. 

• The halo noise from nearby bright stars adds both photon noise and possibly systematic errors in some absorption 
features. For example, we are unable to obtain reliable spectra for bright stars within ~ 0."25 of the WR stars. 

• Photon noise from the background as well as from dark current are significant sources of noise for isolated stars. 

• Spatially varying background, especially emission from Br 7 gas at the Galactic center may lead to systematic 
errors in the equivalent width measurements around Br 7. 



• Read noise from the detector also contributes to the noise, but its contribution is insignificant compared to the 
other sources of error. 
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Some of theses sources of uncertainty affect the early-type and late-type stars differently; for example, the background 
of Br 7 emission may contribute a systematic bias to the intrinsic Br 7 absorption line for young stars, but will not 
affect the measurement of equivalent widths of the Na I doublet. However, the late-type stars will not entirely escape 
this effect because the measurement of the equivalent width of the wavelengths near Br 7 is a parameter in the Bayesian 
evidence. This error is very spatially dependent, as it is the result of a poor estimate of the background at the location 
of the star. The most problema tic regions are therefo re regions where there are strong spatially varying Br 7 emission, 
such as in the Mini-Spiral (e.g. iPaumard et aT1l2004l ). 
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Figure 20. The equivalent width of Br 7 (top) and Na I (bottom) for the sample of stars used for the star planting simulations as a 
function of K' magnitude. 

These complexities result in variations in the sensitivity to early-type and late-type spectra, depending on the 
location of a star. We determine the relative probability that each untyped source is early-type or late-type by running 
a series of Monte Carlo simulations. For each untyped star with brightness K' < 16.0 mag, we simulate and plant 100 
late-type and 100 early-type stars nearby. Each simulated star is planted in the following fashion. 

1. We randomly choose a template spectrum to plant from the catalog of spectra that have already been spectral- 
typed as early or late-type. The template is required to have SNR > 35 to be chosen. The planted young 
stars are also required to have a measured Br 7 equivalent width and excludes all WR stars. These criteria are 
chosen to exclude the spectra from the more massive young stars from the simulations (which may have no Br 
7 absorption). We do not expect the untyped population to include these types of stars as the majority of the 
untyped stars are much fainter. There are 83 late-type and 41 early-type spectra satisfying these requirements. 
Figure [20] shows the equivalent widths of the template spectra used for the simulations. 

2. The spectrum is scaled to the flux corresponding to the magnitude of the untyped sources (the conversion between 
magnitude and flux is empirically calibrated with a power-law fit to the flux determined by Starfinder on the 
OSIRIS cube and the magnitudes determined from NIRC2 imaging). 

3. Photon noise is added to the spectrum. 

4. We then plant the star next to the untyped source at a distance randomly chosen between 4 and 6 pixels in 
radius from the untyped source. The location of the simulated source is also chosen to avoid falling on sources 
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detected in deep imaging. Placing the simulated source- close to the untyped source helps to sample the same 
background location and the same halo noise. See Figure [21] for an example of the locations of simulated sources 
with respect to the untyped sources. 
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Figure 21. Example of a star planting simulation. This example is taken from the field GC East (PA = 285 deg), shown before star 
planting simulations. Labeled with green circles are stars with K 1 < 16.0 and do not have a spectral type assigned by hand (stars not 
in Sample 1). In each simulation, we plant a star next to each of these sources between 4 to 6 pixels from the untyped star (red circles). 
These locations are also chosen so that they will not fall on another source known from imaging or another simulated source. By planting 
stars next to the untyped sources, we can sample the environmental factors that contribute to the incompleteness, such as halo noise from 
being near bright sources. 



UNBINNED FITTING OF RADIAL SURFACE DENSITY PROFILES USING BAYESIAN ANALYSIS 

In order to avoid binning the stellar surface-density profile by radius, we will combine the individual likelihood of 
each source and compute the posterior distribution for the power-law slope. V using Baycs' Theorem: 

P (r|c) = £<™ (E1) 

where D are the data points, P(D\T) is the likelihood, P(T) is the prior distribution of T, and P{D) is the evidence. 
We will assume a flat prior for T. The surface-density profile is used as the likelihood: 



X(x,y,r)oc (v^T^)-r (E2 ) 

where x, y are the RA and DEC projected positional offsets from Sgr A*. We incorporate the individual positional 
measurements and their errors by convolving the density profile by a normalized Gaussian (G(x,y,ax,ay)) centered 
at the measured position with the error as the a. We also normalize the likelihood by integrating over the area of the 
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survey. The likelihood for an individual star is then: 

C(x- y-\T a a ) = J ^{x,y,T)G{x l ,y t \x,y,a Xz ,a yi )dxdy 

V * j ' ^(x,y,T)G{x i ,y i \x,y,(j Xt ,(jy t )dxdydxidx l 
= J y, r)G(xj, y l \x, y, a Xi , a Vi )dxdy 
J y, T)dxdy 

The membership probability of the star (whether it is young or old) can be easily included by raising the likelihood to 
the power of the associated probability (Pe or Pl)', this weighs the likelihood by the star's spectral-type probability. 
In order to incorporate image completeness, we modify weight by the image completeness at the magnitude of the 
given star, I(K'), at the field location: P we i g ht = Ptypel 'I(K'). The total likelihood is a product of all the individual 
likelihoods: 

N 

P(D\T) = C to t a i =Y[£ i (x i ,y i \T,a Xi ,a yi ) p ^ ht (E5) 

i 

The posterior distribution P(T\D) is then sampled using a Markov-Chain Monte Carlo (MCMC) using the Metropolis- 
Hastin gs method. The chains are tested for convergence using the power spectrum method described in lDunklev et akl 
(|2005| ). The Bayesian analysis leads to a natural way of using all available information in determining the stellar 
density power-law and is especially useful for determining the properties of subsamples of stars such as splitting the 
young stars by magnitude. Figure [H] shows the posterior distribution for the completeness-corrected late-type and 
early-type stars, as well as the subpopulation of faint and bright stars. 




Figure 22. The posterior PDFs for the surface-density power-law slope of the late-type (left) and early-type (right) stars. We also 
determine the slopes of the sub-population of bright (K 1 < 14.3, dotted) and faint (K' > 14.3, dashed) sources. The different sub- 
populations appear consistent with each other. The most probable values for T for the different populations are tabulated in Table [3] 
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Table 6 

OSIRIS observations of late-type stars 
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0.09 


11.94 


1.97 


IRS29S 


11.22 
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0.12 


15.16 
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-2.77 


S3-291 


12 


07 


0. 
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11 
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03 


13 


,58 


6.49 


S10-18 


13 


37 


0. 


04 


13 


39 


9.29 


S7-24 


13 


,17 





10 


13 


,57 


6.25 


S10-20 


13 


53 





03 


13 


,78 


10.34 


S9-26 


13 


.56 


0. 


16 


13 


56 


8.80 


S10-21 


13 


66 





11 


13 


,68 


8.84 


S8-21 


13 


.60 


0. 


05 


13 


81 


8.25 


Sll-26 


13 


71 


0. 


1 


13 


83 


11.72 


S10-23 


13 


73 





08 


13 


,77 


9.09 


S6-104 


13 


.94 


0. 


04 


13 


96 


6.69 


S12-19 


13 


91 


0. 


05 


11 


02 


12.00 


S6-106 


13 


96 





04 


13 


99 


6.70 


S9-34 


13 


.99 


0. 


12 


13 


90 


7.37 


Sll-32 


14 


08 





07 


14 


,27 


11.39 


S10-29 


14 


21 





05 


14 


33 


10.96 


S8-36 


14 


.25 


0. 


03 


11 


42 


7.95 


Sll-35 


14 


50 





02 


14 


21 


9.94 


S7-34 


14 


.14 





13 


14 


13 


7.08 


Sll-36 


11 


40 


0. 


05 


11 


50 


10.62 


S10-31 


14 


33 





04 


14 


,11 


10.82 


S12-29 


14 


47 





08 


14 


58 


12.20 


S7-35 


14 


.33 


0. 


12 


11 


26 


7.06 


S9-43 


14 


.39 





09 


14 


29 


7.90 



iC offset 


K 


Epoch 


Na 


(") 
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(A) 


-1.52 


3.74 


2008 


,62 


2.15 


2.54 


3.75 


2007 


12 


5.69 


-1.26 


3.75 


2007 


,88 


5.69 


0.97 


3.96 


2008 


.41 


3.08 


-3.92 


3.93 


2008 


18 


5.52 


-1.48 


3.97 


2007 


,75 


2.89 


-0.36 


4.03 


2008 


.42 


3.96 


1.66 


4.11 


2007 


.90 


5.30 


-1.48 


4.12 


2007 


.96 


1.01 


0.12 


4.22 


2008 


.35 


4.55 


-3.61 


4.25 


2008 


16 


4.93 


-2.74 


4.27 


2008 


.14 


5.19 


-2.23 


4.31 


2007 


.94 


6.73 


-2.43 


4.34 


2007 


.95 


4.16 


-3.61 


4.34 


2008 


,25 


3.65 


-3.25 


4.36 


2007 


.99 


5.22 


-0.26 


4.42 


2008 


.24 


5.83 


-3.16 


4.44 


2008 


.12 


4.85 


-3.75 


4.61 


2008 


.02 


4.09 


0.21 


4.76 


2008 


,75 


4.89 


-1.54 


4.85 


2009 


.04 


2.94 


0.85 


4.88 


2008 


,52 


1.64 


-1.38 


4.88 


2008 


11 


4.46 


-4.79 


5.13 


2008 


.49 


3.81 


-0.94 


5.29 


2008 


.90 


4.88 


-2.74 


5.36 


2008 


17 


5.64 


-2.55 


5.50 


2009 


.31 


6.42 


-2.86 


5.66 


2009 


.03 


2.87 


-3.03 


5.72 


2009 


,28 


2.61 


-3.91 


5.93 


2008 


11 


1.51 


-4.39 


5.96 


2009 


,52 


1.46 


-3.26 


5.95 


2007 


.91 


4.78 


-3.67 


6.13 


2008 


70 


3.47 


-4.68 


6.21 


2008 


15 


4.88 


-4.52 


6.88 


2008 


26 


5.46 


-4.75 


6.93 


2008 


,82 


1.84 


-4.86 


6.98 


2008 


,52 


5.75 


0.90 


6.97 


2008 


.03 


3.68 


-0.20 


6.97 


2008 


15 


4.89 


-1.44 


6.39 


2007 


,89 


4.89 


2.19 


6.17 


2008 


02 


2.70 


1.86 


7.29 


2007 


,98 


5.40 


-0.46 


11.33 


2008 


.49 


5.59 


-2.09 


11.85 


2008 


.01 


5.58 


-0.58 


11.45 


2008 


.10 


1.73 


-6.58 


12.90 


2007 


.95 


6.44 


-0.98 


12.15 


2008 


16 


5.17 


-5.84 


12.00 


2008 


.37 


5.08 


-0.43 


7.53 


2007 


.95 


5.91 


1.57 


7.83 


2007 


.92 


7.44 


0.57 


5.74 


2007 


.88 


4.21 


-5.51 


9.73 


2008 


.20 


4.56 


-4.41 


7.84 


2008 


,08 


4.32 


-3.75 


10.02 


2008 


.23 


4.32 


-3.94 


7.38 


2008 


.19 


3.86 


-0.26 


10.34 


2008 


16 


3.82 


-2.51 


9.15 


2008 


18 


4.45 


-5.09 


10.20 


2007 


.81 


4.50 


0.67 


8.28 


2008 


.31 


4.63 


-1.11 


11.77 


2008 


.12 


3.77 


-4.47 


10.13 


2008 


.49 


5.23 


1.80 


6.92 


2007 


.80 


2.02 


-2.16 


12.19 


2008 


.33 


3.92 


1.38 


6.84 


2008 


.19 


2.95 


-5.17 


9.00 


2008 


.20 


2.79 


0.64 


11.40 


2008 


.46 


3.95 


-0.03 


10.96 


2008 


.32 


3.98 


0.92 


8.01 


2008 


.31 


4.61 


-6.13 


11.68 


2008 


.54 


6.50 


-3.26 


7.79 


2007 


.96 


3.26 


-4.63 


11.58 


2007 


,89 


5.74 


-0.91 


10.85 


2008 


.48 


4.92 


-1.21 


12.26 


2008 


.29 


4.13 


-0.91 


7.12 


2008 


.43 


4.97 


-5.83 


9.82 


2007 


.95 


6.39 



°~Na 

(A) 



N„ 



SNR C Field 



2.72 
0.25 
0.43 

0.08 
1.11 
0.14 
0.29 
2.21 
0.96 
0.61 
0.19 
0.30 
0.28 
0.78 
0.34 
0.17 
0.68 
0.54 
4.33 
0.82 
0.67 
0.33 
1.01 
0.97 
0.86 
0.79 
2.16 
1.17 
0.30 
1.70 
0.41 
1.44 
0.87 
0.21 
0.45 
1.71 
1.75 
0.41 
0.50 
0.74 
0.14 
0.08 
0.39 
0.06 
0.48 
0.49 
0.36 
0.06 
0.59 
1.67 
1.59 
1.00 

0.48 
0.07 
0.55 
0.35 
0.29 
0.96 
0.19 
0.63 
0.88 
0.48 
0.23 

0.22 
0.19 
0.82 
0.48 
0.39 
1.56 
1.65 
0.91 
1.79 
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3 


31 


E2-1 


6 


27 


S2-1 


6 
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4 


10 
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1 


19 


E2-3 


3 


23 


S2-1 


1 


36 


E2-3 


6 


30 
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6 


34 
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4 


29 
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4 


9 


E2-1 


6 


29 
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6 


9 
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6 


37 
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6 


24 
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6 


36 
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1 


25 
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3 


14 
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4 


30 
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6 


22 
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6 


45 
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6 


22 
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6 


51 
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6 


11 
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6 


51 


E2-3 


6 


33 
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6 


48 
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5 


24 


E2-3 


6 


67 


E2-1 


3 


62 


E3-1 


6 


46 


E2-2 


4 


10 
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3 


46 
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6 


62 


E4-1 


3 


39 
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6 


53 


E4-1 


1 


36 


E4-3 


6 


57 
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7 


43 
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6 


15 
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3 


35 


E3-1 


6 


3 1 


E2-1 


6 


17 


E3-3 


6 


40 


E3-3 


1 


23 


E3-2 


6 


33 


E3-3 


4 


39 


E4-1 


6 


58 
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3 


15 
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6 


39 
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6 


19 
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3 


11 
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1 


32 
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3 


11 
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6 


25 
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6 


19 
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1 


19 
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6 


15 
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21 
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•5 


19 
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3 


32 
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Table 6 — Continued 
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S6-112 


14 


30 





09 


14 


37 


6.02 


S12-31 


1 1 


56 





07 


1 1 


16 


12.50 


S12-32 


1 1 


15 





06 


1 1 


56 


11.80 


S9-45 


14 


58 





09 


14 


49 


7.88 


S9-46 


14 


.51 





08 


1 1 


13 


8.06 


S9-47 


11 


18 





03 


1 1 


73 


9.89 


S5-243 


14 


63 





06 


14 


65 


5.52 


S10-41 


1 1 


70 





01 


1 1 


95 


10.09 


S8-47 


14 


.71 





07 


14 


97 


8.79 


S6-122 


1 1 


8 1 





13 


14 


75 


6.71 


S13-32 


15 


06 





06 


1 1 


98 


11.86 


S8-56 


14 


.82 





11 


15 


05 


8.74 


S7-46 


14 


.90 





09 


14 


98 


6.21 


S7-47 


14 


87 





15 


1 1 


8 1 


7.50 


S6-128 


14 


76 





11 


14 


88 


5.71 


S10-54 


15 


03 





04 


15 


29 


10.39 


Sll-49 


15 


39 





01 


15 


2 1 


9.75 


Sll-52 


15 


.14 





07 


15 


26 


11.59 


S10-59 


15 


04 





05 


15 


02 


10.01 


S10-60 


1 1 


92 





01 


15 


18 


10.19 


S10-61 


15 


02 





03 


15 


13 


10.61 


S8-63 


14 


.95 





1 1 


1 1 


.97 


7.59 


S12-47 


15 


12 





06 


15 


2 1 


12.02 


S9-57 


14 


88 





05 


15 


.04 


9.60 


S12-48 


15 


29 





10 


15 


10 


12.40 


S10-65 


15 


10 





06 


15 


10 


9.31 


S12-51 


15 


21 





06 


15 


30 


11.77 


S10-68 


1 1 


91 





05 


14 


92 


9.58 


S10-69 


15 


20 





08 


15 


25 


9.18 


S7-53 


15 


11 





13 


15 


00 


6.91 


S8-65 


15 


.20 





07 


15 


16 


8.90 


Sll-61 


15 


25 





06 


15 


34 


10.79 


S7-58 


15 


.19 





15 


15 


09 


7.12 


S8-66 


15 


.32 





1 1 


15 


10 


8.50 


S7-59 


15 


.20 






15 


.30 


6.02 


S7-61 


15 


.23 





11 


15 


36 


7.86 


S7-62 


15 


22 





1 1 


15 


12 


7.22 


S8-69 


15 


17 





11 


15 


18 


7.57 


Sll-65 


15 


39 





07 


15 


19 


11.29 


S8-72 


15 


27 





06 


15 


50 


8.23 


S7-70 


15 


13 





09 


15 


46 


7.76 


S6-142 


15 


22 





10 


15 


36 


5.46 


S7-72 


15 


21 





10 


15 


19 


6.46 


S9-81 


15 


24 





1 1 


15 


23 


9.18 


S5-247 


15 


3 1 





1 1 


15 


18 


5.66 


S8-76 


15 


17 





13 


15 


18 


7.40 


S9-83 


15 


.29 





10 


15 


22 


7.98 


S9-87 


15 


.33 





1 1 


15 


33 


8.79 


S9-88 


15 


39 





13 


15 


11 


8.52 


S6-148 


15 


45 


0. 


12 


15 


37 


6.20 


S9-93 


15 


■15 





12 


15 


17 


8.82 


S9-95 


15 


.11 





15 


15 


41 


8.97 


S6-152 


15 


.52 





13 


15 


12 


6.57 


S7-83 


15 


38 





11 


15 


33 


7.20 


S9-97 


15 


15 





05 


15 


1 1 


9.27 


S7-85 


15 


.45 





12 


15 


41 


6.61 


S7-88 


15 


.60 





12 


15 


50 


6.69 


Sll-74 


15 


.14 


0. 


06 


15 


25 


11.00 


S10-93 


15 


06 





03 


15 


.07 


10.18 


S10-94 


15 


38 





03 


15 


.45 


9.86 


S13-57 


15 


22 


0. 


13 


15 


37 


13.08 


S12-62 


15 


22 





08 


15 


.19 


12.18 


S8-110 


15 


41 


0. 


21 


15 


.48 


8.63 


S9-102 


15 


42 


0. 


05 


15 


,19 


9.33 


S10-96 


15 


23 





01 


15 


.34 


10.10 


S10-97 


15 


20 


0. 


05 


15 


31 


10.22 


S10-98 


15 


21 


0. 


06 


15 


.29 


9.85 


S8-121 


15 


.45 





08 


15 


46 


6.76 


S12-71 


15 


10 





14 


14 


99 


11.38 


S12-75 


15 


36 





07 


15 


,15 


12.71 


S10-102 


15 


32 





01 


15 


42 


10.79 
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2008.42 


5.53 


0.83 


4 


32 


E2-2 


2008.11 


3.22 


1.01 


4 


25 


E4-2 


2007.85 


4.59 


0.72 


6 


55 


E4-2 


2008.39 


2.85 


1.88 


6 


30 


E3-3 


2008.19 


5.51 


0.09 


6 


26 


E3-3 


2007.91 


6.39 


1.41 


6 


29 


E3-1 


2007.60 


2.61 


1.37 


6 


23 


E2-1 


2008.39 


3.64 


0.55 


6 


37 


E3-1 


2007.90 


2.68 


0.65 


6 


38 


E3-1 


2008.08 


4.28 


2.77 


1 


31 


E3-2 


2007.80 


3.61 


1.27 


1 


22 


E4-3 


2008.10 


2.71 


0.84 


6 


23 


E3-1 


2008.58 


2.66 


1.20 


6 


29 


E3-3 


2008.60 


5.39 


0.59 


6 


35 


E3-2 


2008.91 


4.49 


0.75 


1 


31 


E2-2 


2008.70 


5.65 




2 


15 


E3-1 


2007.70 


2.20 


1.71 


6 


15 


E4-3 


2008.08 


3.99 


0.54 


6 


19 


E4-1 


2007.53 


5.57 


0.86 


6 


22 


E4-2 


2008.47 


2.87 


0.68 


1 


20 


E3-1 


2007.37 


5.44 


1.52 


6 


17 


E4-1 


2008.86 


2.26 


0.37 


3 


20 


E3-3 


2007.94 


5.50 


1.50 


6 


17 


E4-1 


2008.40 


5.15 


0.39 


6 


25 


E3-1 


2008.11 


4.96 


1.62 


6 


15 


E4-1 


2007.54 


3.25 


1.19 


7 


19 


E4-3 


2008.53 


4.92 


1.87 


6 


40 


E4-2 


2007.68 


2.76 




2 


20 


E3-2 


2008.12 


4.20 


2.22 


1 


22 


E4-3 


2007.95 


6.26 


0.77 


6 


20 


E3-2 


2008.59 


2.23 


1.39 


6 


18 


E3-1 


2008.01 


5.19 


1.20 


6 


31 


E4-2 


2008.74 


4.22 


0.85 


6 


34 


E3-2 


2008.25 


2.50 




2 


14 


E3-2 


2009.36 


4.33 


1.66 


6 


17 


E3-3 


2008.51 


2.43 


0.33 


6 


15 


E3-1 


2007.59 


2.51 


0.43 


6 


21 


E3-2 


2007.99 


7.34 


1.66 


1 


12 


E3-2 


2006.99 


5.49 


12.34 


6 


20 


E4-2 


2008.64 


3.32 


0.92 


6 


21 


E3-1 


2008.05 


3.93 


0.94 


6 


23 


E3-2 


2009.25 


3.28 


0.97 


6 


21 


E2-3 


2009.19 


4.98 




1 


21 


E3-2 


2008.10 


3.67 


1.07 


6 


21 


E3-2 


2009.12 


2.38 


2.90 


6 


15 


E2-2 


2007.83 


3.66 


1.07 


6 


21 


E3-3 


2007.88 


5.41 


0.95 


6 


21 


E3-3 


2008.87 


3.48 


1.56 


6 


25 


E3-2 


2008.07 


6.05 


1.13 


6 


25 


E3-3 


2007.61 


5.07 


0.92 


3 


13 


E2-2 


2009.00 


4.67 




2 


18 


E3-2 


2007.66 


3.20 


0.73 


6 


23 


E3-2 


2008.33 


3.96 


0.88 


3 


16 


E3-2 


2008.30 


3.72 


0.87 


6 


32 


E3-2 


2008.60 


4.16 


1.06 


6 


16 


E3-2 


2007.84 


5.92 


1.47 


3 


21 


E3-2 


2007.79 


3.80 


3.01 


1 


11 


E3-2 


2008.42 


3.68 


0.42 


6 


20 


E4-2 


2008.36 


6.28 


0.38 


6 


18 


E4-2 


2008.06 


2.58 


3.03 


6 


12 


E4-2 


2008.77 


3.08 


2.35 


3 


26 


E4-1 


2008.13 


2.68 


0.46 


6 


21 


E4-2 


2008.57 


3.71 


0.33 


6 


36 


E3-2 


2008.69 


2.64 




2 


18 


E3-2 


2008.14 


3.73 


1.58 


6 


21 


E4-3 


2007.48 


3.25 


2.15 


5 


17 


E4-2 


2007.28 


3.92 


2.48 


5 


17 


E4-2 


2007.23 


2.95 


0.92 


6 


21 


E3-3 


2007.34 


4.64 


2.90 


6 


17 


E4-3 


2008.31 


3.34 


1.99 


6 


21 


E4-1 


2008.26 


2.85 


0.99 


6 


11 


E4-1 



-0.59 6.05 

-3.43 12.97 

-2.50 12.06 

-5.43 9.57 

-5.41 9.70 

-0.26 9.89 

1.45 5.70 

0.65 10.11 

0.65 8.81 

-1.74 6.93 

-5.97 13.28 

0.08 8.74 

-4.25 7.52 

-2.75 7.99 

-2.53 6.25 

0.59 10.40 

-5.82 11.35 

-1.20 11.66 

-2.31 10.27 

0.93 10.23 

-0.31 10.61 

-3.37 8.30 

-1.60 12.12 

-0.49 9.61 

-1.63 12.51 

-5.24 10.68 

-3.11 12.18 

-3.04 10.05 

-4.70 10.32 

-2.45 7.33 

0.71 8.93 

-2.44 11.06 

-1.37 7.25 

-0.82 8.54 

-4.84 7.73 

-0.10 7.86 

-1.85 7.45 

-3.25 8.23 

-4.00 11.97 

1.21 8.32 

-1.13 7.84 

-3.56 6.52 

-3.10 7.17 

-2.63 9.55 

-1.65 5.89 

-3.82 8.32 

-4.85 9.34 

-2.37 9.11 

-4.06 9.44 

-2.38 6.64 

-3.54 9.50 

-3.15 9.51 

-1.97 6.86 

-2.81 7.73 

-2.39 9.57 

-2.81 7.19 

-2.44 7.12 

-2.37 11.25 

-2.75 10.54 

-3.62 10.50 

-0.06 13.08 

-3.31 12.63 

-1.58 8.77 

-1.13 9.40 

-4.33 10.97 

-3.91 10.92 

-3.74 10.53 

-4.92 8.33 

-6.15 12.90 

-1.64 12.79 

-1.17 10.85 



Table 6 — Continued 



Name 


K' 




K AA RA offsct DEC offsct 


R 


Epoch 


Na 


o-Na N obB b SNR C Field 








(") (") 


(") 




(A) 


(A) 



Corrected for differential extinction to an Ak s — 2.7 from lSchodel et all fl2010l) . 

Each observation has an integration time of 900 s. 

Signal to noise per pixel calculated between 2.212 to 2.218 fj,m. 

Identified as early-type in lDo et a l. (2009a), but now as late-type. See Appcndix^^f 01 " rnore details 



Table 7 

OSIRIS observations of early-type stars 



Name 


K' 




~K' a 
AA 


RA offset 


DEC offset 


R 


Epoch 


Br 7 




N ob3 h 


SNR C 


Field 








(") 


(") 


(") 




(A) 


(A) 








SO-2 


14.13 


0.07 


14.39 


-0.01 


0.17 


0.17 


2007.89 


3.79 


0.29 


2 1 


74 


C 


S0-1 


14.74 


0.07 


14.96 


0.04 


-0.26 


0.26 


2006.30 


1.13 


0.51 


7 


31 


C 


SO-3 


14.54 


0.03 


14.82 


0.34 


0.12 


0.36 


2008.39 


2.37 


0.48 


7 


26 


c 


SO-5 


15.06 


0.12 


15.27 


0.18 


-0.37 


0.41 


2007.99 


7.80 


1.10 


9 


15 


c 


SO- 11 


15.13 


0.04 


15.41 


0.50 


-0.06 


0.51 


2008.08 


10.20 


7.60 


9 


21 


c 


SO-4 


14.31 


0.14 


14.51 


0.42 


-0.31 


0.52 


2007.62 


3.10 


3.00 


9 


37 


c 


SO-9 


14.25 


0.08 


14.42 


0.21 


-0.59 


0.62 


2008.03 


5.40 


2.00 


7 


38 


c 


SO-31 


15.09 


0.10 


15.31 


0.55 


0.45 


0.71 


2007.97 


6.01 


1.54 


6 


22 


c 


S0-14 


13.72 


0.14 


13.88 


-0.76 


-0.28 


0.81 


2008.24 


4.52 


0.19 


1 


17 


c 


Sl-3 


12.09 


0.05 


12.22 


0.36 


0.88 


0.95 


2007.36 


1.31 


0.22 


3 


90 


c 


SO-15 


13.60 


0.09 


13.69 


-0.96 


0.21 


0.98 


2007.77 


1.38 


0.16 


7 


50 


G 


Sl-2 


14.83 


0.17 


14.88 


0.05 


-1.02 


1.02 


2007.90 


3.29 


1.13 


3 


25 


c 


Sl-1 


13.07 


0.03 


13.22 


1.03 


0.03 


1.03 


2007.79 






5 


57 


c 


Sl-4 


12.55 


0.12 


12.67 


0.86 


-0.66 


1.08 


2007.63 






5 


52 


c 


Sl-8 


14.22 


0.13 


14.29 


-0.61 


-0.90 


1.09 


2007.37 


2.68 


0.25 


4 


28 


c 


Sl-33 


14.92 


0.04 


15.05 


-1.24 


-0.01 


1.25 


2008.07 


4.82 


0.24 


1 


13 


c 


Sl-12 


13.57 


0.15 


13.63 


-0.78 


-1.02 


1.28 


2007.72 


2.15 


0.86 


3 


20 


G 


Sl-14 


12.73 


0.06 


12.88 


-1.33 


-0.36 


1.38 


2008.16 


1.22 


0.20 


2 


31 


c 


Sl-22 


12.57 


0.05 


12.65 


-1.59 


-0.51 


1.67 


2007.80 






3 


61 


c 


Sl-19 


13.62 


0.11 


13.58 


0.42 


-1.63 


1.68 


2008.11 


7.00 


1.00 


1 


38 


s 


Sl-18 


14.84 


0.06 


14.75 


-0.78 


1.51 


1.69 


2008.47 


1.81 


1.00 


2 


2 1 


N 


Sl-24 


11.49 


0.12 


11.47 


0.73 


-1.64 


1.79 


2007.87 


2.84 




1 


71 


SE 


S2-7 


14.21 


0.13 


14.11 


0.94 


1.85 


2.07 


2008.01 






5 


46 


N 


IRS16CC 


10.68 


0.07 


10.75 


1.99 


0.59 


2.07 


2007.92 










NE 


S2-5 


13.29 


0.09 


13.32 


1.93 


-0.79 


2.09 


2008.23 






7 


58 


E 


S2-4 


12.09 


0.12 


12.09 


1.50 


-1.46 


2.10 


2008.27 


1.94 


0.19 


11 


149 


SE 


S2-6 


12.04 


0.14 


12.02 


1.64 


-1.34 


2.12 


2007.92 


2.04 


0.16 


7 


103 


SE 


S2-40 


15.37 


0.04 


15.42 


1.73 


1.28 


2.15 


2008.16 






5 


16 


NE 


IRS33N 


11.29 


0.12 


11.23 


-0.03 


-2.23 


2.23 


2008.08 


3.10 


0.40 


10 


115 


S 


S2-50 


15.42 


0.13 


15.38 


1.70 


-1.51 


2.27 


2008.49 


8.00 


3.00 


11 


16 


SE 


S2-17 


10.74 


0.08 


10.64 


1.32 


-1.88 


2.30 


2007.92 


3.19 


0.07 


11 


125 


SE 


S2-22 


12.86 


0.01 


12.92 


2.30 


-0.22 


2.31 


2007.97 






10 


96 


E 


S2-21 


13.36 


0.09 


13.34 


-1.64 


-1.66 


2.33 


2008.06 


1.20 


0.70 


8 


73 


S 


S2-19 


12.60 


0.04 


12.50 


0.40 


2.31 


2.34 


2007.83 






5 


60 


N 


S2-58 


14.08 


0.10 


14.07 


2.15 


-1.14 


2.43 


2008.26 


4.70 


3.00 


7 


41 


E 


S2-74 


13.16 


0.03 


13.07 


0.14 


2.78 


2.78 


2007.64 






5 


60 


N 


S2-76 


15.07 


0.06 


14.99 


-0.23 


2.81 


2.82 


2008.41 


1.18 


1.00 


5 


28 


N 


S2-29 


15.32 


0.06 


15.25 


1.95 


-2.15 


2.90 


2008.09 






11 


19 


SE 


S3-2 


12.02 


0.04 


12.13 


3.08 


0.55 


3.12 


2007.84 


0.49 


0.04 


3 


17 


NE 


JO-O 




u.uo 


J-O. i-O 


o.uo 


-u.uo 


O. lO 


9007 79 
zuu / . / z 


1 on 


9 nn 
z.uu 


1 n 


21 




S3- 17 


13.55 


0.06 


13.46 


-1.41 


2.85 


3.18 


2007.84 






6 


87 


NW 


S3-96 


14.31 


0.07 


14.21 


-3.13 


-0.64 


3.20 


2007.27 






1 


29 


sw 


S3-30 


12.39 


0.06 


12.30 


1.66 


-2.94 


3.38 


2007.90 


1.63 


0.44 


8 


56 


SE 


IRS13E1 


10.62 


0.03 


10.60 


-2.97 


-1.65 


3.40 


2007.80 


3.98 


0.31 


1 


113 


SW 


S3-190 


13.96 


0.13 


13.72 


-3.18 


1.42 


3.49 


2008.00 


1.20 


1.10 


11 


55 


vv 


S3-10 


12.10 


0.03 


12.21 


3.34 


-1.12 


3.52 


2008.17 






10 


81 


E 


S3- 12 


11.37 


0.37 


11.28 


2.37 


-2.73 


3.61 


2003.37 






11 


111 


SE 


S4-12 


14.67 


0.27 


14.57 


-2.84 


2.84 


4.02 


2008.08 






6 


50 


NW 


S4-71 


12.34 


0.09 


12.28 


0.77 


-4.08 


4.15 


2008.06 






6 


47 


S2-1 


S4-169 


13.49 


0.02 


13.70 


4.42 


0.27 


4.43 


2008.06 


5.72 


1.54 


3 


30 


E2-1 


S4-196 


14.37 


0.05 


14.24 


2.24 


-3.93 


4.52 


2008.65 






6 


32 


S2-1 


S4-287 


13.68 


0.02 


13.82 


0.13 


-4.77 


4.77 


2008.21 


1.21 


0.86 


5 


12 


S2-1 


IRS1W 


10.93 


0.38 


11.09 


5.26 


0.61 


5.29 


2007.79 


1.47 


0.30 


6 


69 


E2-1 


S5-237 


13.21 


0.10 


13.33 


5.50 


1.00 


5.59 


2007.78 


1.72 


0.89 


6 


26 


E2-1 


S5-233 


12.32 


0.03 


12.46 


5.60 


0.67 


5.64 


2008.14 






6 


31 


E2-1 


S5-183 


11.56 


0.03 


11.75 


4.60 


-3.44 


5.74 


2008.18 


2.49 


0.32 


6 


57 


E2-3 


S5-231 


11.99 


0.05 


12.10 


5.81 


0.09 


5.82 


2007.93 






6 


57 


E2-1 


S6-81 


11.02 


0.06 


11.05 


6.36 


0.26 


6.37 


2007.68 


4.02 


0.41 


6 


111 


E2-1 


S7-30 


13.88 


0.12 


13.83 


6.47 


-2.69 


7.01 


2007.87 


3.57 


0.87 


3 


37 


E3-2 


S7-31 


13.98 


0.09 


13.98 


7.27 


-0.29 


7.28 


2008.18 






1 


49 


E3-1 


S7-36 


14.38 


0.10 


14.46 


6.36 


-4.42 


7.75 


2008.49 


0.82 


0.48 


6 


38 


E3-3 


S8-70 


15.33 


0.16 


15.43 


8.05 


-0.78 


8.09 


2008.95 


2.98 




2 


16 


E3-1 



Table 7 — Continued 



Name 


K' 




K' a 
AA 


RA offset 


DEC offset 


R 


Epoch 


Br 7 


<7Br 7 




SNR C 


Field 








(") 


(") 


(") 




(A) 


(A) 








S9-15 


12.97 


0.07 


13.23 


9.02 


0.42 


9.04 


2007.72 






6 


62 


E3-1 


S9-11 


12.72 


0.03 


12.98 


9.44 


1.04 


9.50 


2008.07 






6 


82 


E3-1 


S10-50 


14.74 


0.03 


14.76 


9.59 


-3.17 


10.10 


2008.50 


2.22 


0.72 


1 


31 


E4-2 


S10-32 


14.36 


0.08 


14.35 


10.20 


-1.70 


10.34 


2007.66 


2.35 


1.28 


5 


58 


E4-2 


S10-34 


14.46 


0.09 


14.36 


8.88 


-5.63 


10.51 


2008.14 


1.02 


1.62 


1 


25 


E4-3 



Corrected for differential extinction to an Ak s — 2.7 from Schodcl ct al. (2010) 

Each observation has an integration time of 900 s. 

Signal to noise per pixel calculated between 2.212 to 2.218 /Am. 



Table 8 

Wolf-Rayet stars in the GCOWS field 



Name K' a K , K' AA a RA offset DEC offset R Epoch Field 

(") (") (") 



IRS16C 


9.94 


0.02 


10.10 


1.07 


0.54 


1 


19 


2007.78 


C 


IRS16NW 


10.13 


0.03 


10.14 


0.07 


1.22 


1 


22 


2007.81 


N 


IRS16SW 


10.08 


0.08 


10.15 


1.10 


-0.95 


1, 


15 


2007.87 


C 


IRS16SW-E 


11.13 


0.01 


11.09 


1.88 


-1.12 


2 


19 


2007.99 


E 


S2-16 


11.99 


0.07 


11.69 


-1.05 


2.06 


2 


31 


2007.80 


NW 


IRS16NE 


9.14 


0.05 


9.27 


2.89 


0.98 


3 


05 


2007.84 


NE 


S3-5 


12.00 


0.06 


12.01 


2.95 


-1.16 


3 


17 


2008.18 


E 


IRS33E 


10.20 


0.08 


10.14 


0.69 


-3.13 


3 


21 


2008.12 


S2-1 


IRS13E4 


11.73 


0.06 


11.65 


-3.23 


-1.41 


3 


52 


2008.01 


SW 


IRS13E2 


10.69 


0.04 


10.62 


-3.19 


-1.73 


3 


63 


2007.84 


SW 


S6-82 


13.51 


0.13 


13.48 


6.72 


-0.47 


6 


73 


2008.08 


E2-2 


S9-1 


12.65 


0.01 


12.87 


9.45 


0.28 


9 


45 


2008.47 


E3-1 



Corrected for differential extinction to an -Ak — 2.7 from lSchodcl et al.l ( 1201011 



31 



Table 9 

Completeness Simulation Results 



£1 TTl O 
IN dUlL. 


K' 


CT A" 


K' 

AA 


A RA 


A Dec 


ln(BF) a 


SNR 


Na F,W 


Na PW Prr 

lid J_j VV J_j1 1 


Rr -v F,W 

j_>i y j_j vv 


Rr 'V PW Prr 

JJI V 11/ V V HiL 1 


P(01d) b 




Field 










( ) 








/ X \ 

(, A ) 




( A > 


( A J 








SO-35 


15.30 





,12 


15.48 


0.04 


0.90 


5.86 


17 


5.44 


1.44 


2.86 


2.43 


1.00 


0.00 


c 


Sl-27 


15.35 





.06 


15.43 


-1.03 


0.20 


-3.31 


30 


-0.13 


0.74 


2.36 


0.45 


0.31 


0.69 


c 


Sl-29 


15.28 





01 


15.45 


1.08 


0.16 


4.20 


8 


6.35 


1.85 


4.11 


1.21 


1.00 


0.00 


c 


Sl-32 


15.19 





12 


15.30 


-0.98 


-0.65 


9.66 


20 


5.30 


1.10 


0.83 


0.36 


1.00 


0.00 


c 


Sl-50 


15.36 





.06 


15.40 


1.49 


0.66 


3.09 


11 


4.41 


2.01 


1.10 


0.58 


1.00 


0.00 


c 


Sl-51 


14.97 





07 


15.07 


-1.66 


-0.17 















0.50 


0.50 


c 


Sl-53 


15.20 





.03 


15.28 


1.67 


-0.11 


2.76 


12 


4.14 


1.98 


1.72 


0.55 


1.00 


0.00 


E 


Sl-45 


15.21 





.03 


15.14 


-1.27 


1.09 


1.75 


17 


2.46 


1.18 


2.57 


0.65 


0.96 


0.04 


c 


Sl-54 


15.41 





.03 


15.38 


-1.51 


0.76 


2.16 


6 


7.94 


3.72 


1.30 


1.73 


0.99 


0.01 


w 


Sl-52 


15.13 





19 


15.07 


0.02 


1.69 


1.20 


19 


1.76 


1.02 


1.30 


0.36 


1.00 


0.00 


N 


Sl-49 


14.37 





.06 


14.46 


-1.70 


0.16 




21 










0.50 


0.50 


w 


Sl-55 


15.41 





.05 


15.46 


1.59 


0.63 




8 










0.50 


0.50 


c 


Sl-56 


15.52 





.08 


15.36 


-1.10 


1.33 


4.21 


1 1 


5.34 


1.95 


1.09 


0.33 


1.00 


0.00 


c 


Sl-58 


15.51 





01 


15.41 


-1.48 


0.98 


0.33 


5 


13.32 


7.14 


8.85 


3.95 


0.87 


0.13 


w 


Sl-59 


15.40 





02 


15.34 


0.02 


1.78 


5.51 


15 


4.28 


1.27 


1.41 


0.58 


1.00 


0.00 


N 


Sl-62 


15.35 





01 


15.27 


0.46 


1.79 


2.91 


12 


3.56 


1.56 


1.79 


0.82 


0.98 


0.02 


N 


Sl-63 


15.53 





,25 


15.26 


-1.21 


1.41 


1.61 


1 


4.61 


3.24 


-6.08 


2.71 


0.99 


0.01 


c 


Sl-64 


15.39 





.06 


15.31 


0.66 


1.82 


5.12 


21 


3.89 


1.19 


1.68 


0.42 


1.00 


0.00 


N 


S2-36 


13.38 





06 


13.45 


2.00 


0.43 


-0.15 


24 


1.24 


0.62 


2.25 


0.42 


1.00 


0.00 


E 


S2-43 


15.48 





01 


15.47 


-1.83 


-1.14 


0.96 


11 


2.17 


2.07 


1.22 


1.09 


0.98 


0.02 


sw 


S2-42 


15.35 





11 


15.25 


0.49 


2.14 


4.99 


22 


3.79 


1.17 


1.11 


0.86 


1.00 


0.00 


N 


S2-195 


15.67 





,05 


15.40 


-2.02 


1.46 


-2.16 


7 


-8.71 


3.88 


-4.69 


1.51 


0.12 


0.88 


W 


S2-60 


15.30 





06 


14.90 


-1.44 


2.06 


3.71 


21 


3.89 


1.39 


2.44 


0.69 


1.00 


0.00 


NW 


S2-63 


15.31 





,08 


15.22 


-0.65 


2.48 


1.44 


11 


5.01 


4.34 


1.48 


0.55 


0.99 


0.01 


N 


S2-208 


15.62 





,14 


15.34 


-2.03 


1.59 


-0.66 


8 


0.89 


2.21 


3.21 


1.30 


0.66 


0.34 


NW 


S2-68 


15.34 





10 


14.92 


-1.79 


1.96 


1.85 


15 


4.11 


2.09 


2.96 


0.75 


0.97 


0.03 


NW 


S2-261 


15.51 





02 


15.37 


-2.60 


1.01 


-2.02 


12 


2.72 


1.88 


7.91 


1.51 


0.13 


0.87 


W 


S2-268 


15.42 





03 


15.31 


-2.77 


0.49 




8 










0.61 


0.39 


w 


S2-277 


14.68 





,08 


14.64 


-2.44 


-1.43 


-1.55 


15 


0.06 


1.17 


0.93 


0.59 


0.69 


0.31 


SW 


S2-82 


15.23 





03 


15.30 


2.85 


0.05 


6.43 


28 


4.59 


1.19 


2.24 


0.73 


1.00 


0.00 


E 


S2-79 


15.06 





02 


15.13 


2.86 


-0.13 


1.37 


31 


2.73 


1.12 


3.77 


0.34 


0.99 


0.01 


E 


S2-30 


15.08 





03 


15.16 


2.92 


-0.05 


-1.47 


20 


1.67 


1.08 


4.42 


0.42 


0.81 


0.19 


E 


S2-306 


15.42 





,19 


15.48 


-0.48 


-2.88 


-2.46 


21 


-2.32 


1.73 


1.21 


0.38 


0.69 


0.31 


S 


S2-84 


15.28 





,08 


15.24 


1.66 


-2.47 


12.72 


32 


5.43 


0.92 


1.44 


0.38 


1.00 


0.00 


SE 


S2-86 


15.35 





26 


15.32 


2.67 


-1.37 


0.40 


15 


1.71 


1.30 


1.88 


1.53 


0.97 


0.03 


E 


S3-50 


15.51 





,05 


15.45 


-1.88 


-2.35 


-2.68 


20 


-0.20 


1.21 


3.00 


1.08 


0.53 


0.47 


s 


S3-43 


15.33 





06 


15.36 


-0.17 


-3.01 


3.03 


1 1 


3.92 


1.51 


2.91 


0.85 


1.00 


0.00 


S2-1 


S3-43 


15.33 





06 


15.36 


-0.17 


-3.01 


3.53 


8 


4.85 


1.61 


3.98 


1.47 


1.00 


0.00 


s 


S3-51 


15.11 





03 


15.02 


-0.15 


3.03 


3.15 


14 


3.05 


1.14 


-1.23 


0.56 


1.00 


0.00 


N 


S3-4 


14.74 





,09 


14.80 


3.07 


-0.49 


-1.21 


20 


0.41 


1.05 


1.03 


0.31 


0.68 


0.32 


E 


S3-86 


14.86 





,10 


14.80 


2.11 


-2.31 


4.83 


10 


2.60 


0.68 


0.68 


0.22 


1.00 


0.00 


SE 


S3-92 


15.32 





01 


15.25 


-1.33 


2.88 


-0.27 


18 


1.45 


1.22 


2.89 


0.77 


0.55 


0.45 


NW 


S3-21 


15.38 





03 


15.48 


3.20 


-0.22 


4.26 


22 


3.46 


1.20 


1.24 


0.33 


0.99 


0.01 


E 


S3-23 


15.36 





04 


15.35 


2.90 


-1.38 


5.57 


13 


5.07 


1.45 


2.32 


0.65 


0.99 


0.01 


E 


S3-104 


15.46 





10 


15.35 


-3.16 


-0.66 


-4.60 


13 


1.30 


1.59 


6.29 


0.59 


0.52 


0.48 


SW 


S3-125 


15.42 





,11 


15.37 


-3.07 


-1.22 


-0.52 


<) 


1.02 


5.36 


5.00 


2.41 


0.72 


0.28 


SW 


S3-136 


14.42 





06 


14.38 


-3.02 


-1.41 




4 










0.64 


0.36 


SW 


S3-146 


14.07 





,09 


14.20 


-0.29 


-3.35 


4.35 


12 


7.68 


2.47 


1.86 


0.50 


1.00 


0.00 


S2-1 


S3-403 


14.65 





,05 


14.76 


-0.11 


-3.37 


0.04 


15 


3.04 


1.02 


6.26 


0.85 


0.88 


0.12 


S2-1 


S3-167 


15.37 





,07 


15.16 


-3.07 


1.42 


1.96 


11 


4.64 


2.25 


3.09 


1.32 


0.98 


0.02 


W 


S3-159 


15.17 





,05 


15.10 


-1.28 


3.13 


-2.07 


14 


-0.58 


1.75 


3.01 


0.48 


0.31 


0.69 


NW 


S3-156 


12.99 





,11 


12.97 


0.42 


-3.36 


0.61 


25 


1.75 


0.74 


-1.46 


0.20 


0.89 


0.11 


S2-1 


S3-162 


14.52 





,08 


14.62 


-0.08 


-3.41 


-0.41 


14 


3.25 


1.14 


6.72 


0.92 


0.80 


0.20 


S2-1 


S3-169 


15.46 





,04 


15.19 


-3.39 


0.38 


2.46 


16 


2.98 


1.24 


2.75 


0.70 


0.97 


0.03 


W 


S3-31 


15.35 





10 


15.48 


3.39 


0.38 


-4.14 


8 


-0.15 


1.84 


6.79 


1.28 


0.20 


0.80 


NE 


S3-33 


15.30 





06 


15.45 


3.33 


-0.84 


4.19 


17 


4.54 


1.46 


2.96 


0.67 


1.00 


0.00 


E 


S3-172 


15.35 





06 


15.48 


3.39 


-0.59 


8.04 


38 


4.44 


1.01 


1.43 


0.63 


1.00 


0.00 


E 


S3-192 


15.11 





,07 


15.10 


0.48 


-3.47 


4.44 


9 


5.39 


1.81 


-1.08 


0.79 


1.00 


0.00 


S2-1 


S3-200 


15.54 





,09 


15.43 


-3.31 


-1.17 




1 










0.65 


0.35 


SW 


S3-216 


14.95 





06 


14.66 


-3.13 


1.68 




8 










0.66 


0.34 


W 


S3-229 


15.56 





11 


15.47 


1.56 


-3.23 




7 










0.66 


0.34 


S2-1 


S3-223 


14.90 





04 


14.81 


1.27 


-3.37 


-1.26 


15 


4.15 


1.60 


7.86 


1.12 


0.80 


0.20 


S2-1 


S3-227 


15.27 





13 


15.39 


0.02 


-3.60 


2.91 


9 


4.15 


1.85 


-0.68 


0.42 


0.99 


0.01 


S2-1 


S3-279 


14.74 





13 


14.59 


-3.42 


-1.51 




1 










0.66 


0.34 


SW 


S3-338 


15.75 





02 


15.48 


-3.82 


-0.44 


0.91 


10 


2.93 


1.83 


3.86 


2.18 


0.95 


0.05 


SW 


S3-348 


15.21 





,05 


15.40 


3.90 


0.23 




8 










0.43 


0.57 


E2-1 


S3-348 


15.21 





,05 


15.40 


3.90 


0.23 


-4.31 


8 


-2.73 


1.50 


3.12 


0.80 


0.43 


0.57 


E2-2 


S4-45 


15.20 





20 


15.40 


4.10 


-0.42 


-0.62 


14 


3.91 


1.60 


5.64 


0.49 


0.96 


0.04 


E2-2 


S4-86 


15.22 





,07 


15.40 


-0.67 


-4.12 


0.18 


16 


1.75 


1.58 


2.67 


1.19 


0.99 


0.01 


S2-1 


S4-170 


14.16 





,11 


14.33 


-0.57 


-4.41 


11.28 


25 


5.11 


0.88 


3.16 


1.73 


1.00 


0.00 


S2-1 


S4-280 


15.52 





,04 


15.48 


1.21 


-4.60 


3.53 


4 


5.33 


1.78 


3.40 


0.59 


1.00 


0.00 


S2-1 


S4-315 


15.00 





,07 


15.19 


4.81 


0.84 


3.26 


11 


3.42 


1.22 


2.58 


3.14 


1.00 


0.00 


E2-1 


S5-60 


15.11 





,05 


15.31 


5.21 


0.13 


-1.02 


5 


-0.02 


4.61 


-7.00 


1.57 


0.76 


0.24 


E2-1 


S5-98 


15.16 





,04 


15.29 


1.05 


-5.24 




13 










0.73 


0.27 


S2-1 
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Table 9 — Continued 



Name 


K' 


cr K , 


K' A , 

AA 


A RA 

\ 1 


A Dec 

\ I 


ln(BF) a 


SNR 


Na EW 


Na EW Err 

(A) 


Br 7 EW 
(A) 


Br 7 EW Err 
(A) 


P(01d) b 


P(Yng) b 


Field 


S5-240 


13.92 





07 


14.05 


5.35 


0.10 


1.53 


14 


2.27 


0.99 


-3.37 


1.14 


1.00 


0.00 


E2-1 


S5-94 


15.46 





.07 


15.49 


4.92 


2.11 


3.15 


17 


2.89 


1.13 


-0.39 


0.30 


1.00 


0.00 


E2-1 


S5-256 


15.38 


o 


.14 


15.49 


5.63 


0.87 




1 










0.74 


0.26 


E2-1 


S5-241 


14.00 


o 


.04 


14.09 


5.85 


0.41 


-3.26 


9 


-1.31 


2.03 


4.49 


0.89 


0.51 


0.49 


E2-1 


S6-137 


15.26 


o 


.04 


15.35 


6.01 


0.46 


-1.59 


7 


-2.48 


2.60 


1.42 


1.91 


0.53 


0.47 


E2-1 


S6-129 


15.00 


o 


08 


15.07 


6.10 


-0.33 


-2.17 


10 


-1.92 


2.16 


2.82 


0.96 


1.00 


0.00 


E2-1 


S6-129 


15.00 


o 


08 


15.07 


6.10 


-0.33 


3.78 


22 


2.94 


1.05 


0.46 


0.43 


1.00 


0.00 


E2-2 


S6-131 


14.95 


o 


10 


15.00 


6.18 


0.43 


-1.52 


9 


0.11 


2.09 


3.94 


1.52 


0.88 


0.12 


E2-1 


S6-111 


14.27 


o 


.02 


14.33 


6.30 


0.83 




29 










0.76 


0.24 


E2-1 


S6-138 


15.30 


o 


11 


15.29 


5.93 


-2.32 


3.18 


25 


2.73 


0.84 


3.27 


0.80 


1.00 


0.00 


E2-2 


S6-132 


14.86 


o 


11 


14.86 


6.49 


-0.06 


0.44 


12 


1.63 


1.30 


0.28 


2.02 


0.97 


0.03 


E2-1 


S6-145 


15.31 


o 


08 


15.40 


5.92 


-3.23 


8.72 


16 


5.74 


1.27 


1.19 


0.84 


1.00 


0.00 


E2-3 


S6-151 


15.23 


o 


11 


15.25 


6.81 


0.26 




6 










0.77 


0.23 


E2-1 


S7-81 


15.32 


o 


09 


15.32 


6.34 


-3.19 


0.97 


11 


4.16 


1.78 


4.53 


0.57 


0.99 


0.01 


E3-3 


S7-77 


15.35 


o 


10 


15.37 


7.21 


0.28 


2.04 


15 


2.76 


1.32 


2.35 


0.87 


1.00 


0.00 


E3-1 


S7-65 


15.22 





.08 


15.31 


7.14 


1.29 


1.72 


7 


5.65 


2.46 


4.18 


1.26 


1.00 


0.00 


E2-1 


S7-99 


15.33 


o 


10 


15.36 


6.63 


-3.19 


0.70 


26 


2.69 


1.51 


4.14 


1.32 


0.98 


0.02 


E3-3 


S7-110 


14.98 


o 


20 


15.02 


7.44 


0.33 


4.43 


38 


3.38 


0.73 


5.53 


0.95 


1.00 


0.00 


E3-1 


S7-66 


15.41 





10 


15.47 


6.73 


-3.99 


2.44 


28 


2.46 


1.02 


2.14 


0.32 


1.00 


0.00 


E3-3 


S8-86 


15.31 





10 


15.32 


7.49 


-3.42 


1.11 


21 


2.32 


2.27 


1.68 


0.38 


0.98 


0.02 


E3-3 


S9-91 


15.47 





.09 


15.44 


8.30 


-4.82 


1.70 


18 


2.73 


1.33 


-1.47 


0.42 


1.00 


0.00 


E3-3 


S9-72 


15.40 





.05 


15.32 


8.21 


-5.52 


-1.20 


18 


-1.09 


2.14 


1.90 


0.55 


0.82 


0.18 


E3-3 


S10-46 


14.68 





.12 


14.73 


9.03 


-4.45 


12.64 


50 


5.04 


0.85 


0.93 


0.24 


1.00 


0.00 


E4-3 


S10-319 


15.30 





03 


15.30 


10.12 


-2.79 




1 










0.83 


0.17 


E4-2 


S10-80 


15.24 





.12 


15.29 


9.50 


-5.06 


1.94 


25 


1.95 


0.83 


1.41 


0.44 


1.00 


0.00 


E4-3 


S10-72 


15.32 





.09 


15.41 


9.75 


-4.79 


-2.78 


26 


-0.49 


1.03 


1.66 


0.57 


0.24 


0.76 


E4-3 


Sll-47 


15.51 





11 


15.25 


10.00 


-5.94 




7 










0.84 


0.16 


E4-3 


Sll-60 


15.30 





.08 


15.42 


11.88 


-0.04 


-2.26 


17 


0.91 


1.51 


5.28 


1.27 


0.59 


0.41 


E4-1 


S12-66 


15.26 





06 


15.37 


11.89 


-2.55 




1 










0.85 


0.15 


E4-2 


S12-29 


14.47 





.08 


14.58 


12.20 


-1.22 


3.86 


24 


3.74 


1.18 


-5.34 


1.46 


1.00 


0.00 


E4-1 


S12-211 


15.59 





.06 


15.48 


11.21 


-5.81 


0.92 


11 


3.10 


4.46 


1.80 


0.59 


1.00 


0.00 


E4-3 



Blank entries arc for stars for which wc arc unable to measure a Baycs factor, so the probabilities will be the prior. 

Probability used in the analysis for this paper. This probability is calculated from Equations 1 10 1 and 1 1 1 1 which includes our relative sensitivity to the two types of 
stars as well as a prior on the relative radial surface-density profiles. 



